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In  the  Stanford  9.1  cm  synoptic  charts  during  the  same  period,  when  clear  X-ray 
data  of  coronal  holes  were  available  from  Skylab. 

Chapter  2  summarizes  our  efforts  to  Identify  Individual  high  latitude corona] 
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coronal  holes  were  Identified  from  the  Kitt  Peak  10830  fi  synoptic  maps,  while  the 
high  speed  solar  «rlnd  streams  were  Identified  from  the  interplanetary  scintilla¬ 
tion  measurements  made  at  UCSD  and  at  the  Toyokawa  Observatory  In  Japan. 

Chapter  3  describes  our  1977  coronal  hole  observations  with  the  Haystack 
(1.9A  and  3.81  cm)  and  Areclbo  (11.5  and  21.4  cm)  radio  telescopes./  The  Areclbo 
observations  were  of  special  Importance  because  they  demonstrated  t^t  In  the 
radio  domain  coronal  holes  become  most  clearly  visible  In  the  10-10|0  cm  range, 
provided  the  radio  telescope  has  sufficient  angular  resolution.  Tl^se  observa¬ 
tions  produced  a  radio  spectrum  of  a  coronal  hole. 

Chapter  4  describes  Che  International  effort  we  organized  In  the  summer  of 
1978  to  obtain  solar  radio  maps  over  a  wide  range  of  wavelengths.  These  radio 
maps  were  made  on  the  same  days  at  a  number  of  radio  observatories  around  the 
world.  ')  During  the  first  observing  period,  July  6-10,  1978  there  were  no  coronal 
holes  {/resent.  During  the  second  period,  July  15-18,  1978,  there  were  two  siz¬ 
able  coronal  holes,  one  In  the  northern  hemisphere  and  the  other  In  the  southern 
hemls]^here.  Unfortunately  there  were 'also  several  active  regions,  which  made  It 
difficult  to  establish  the  quiet  sun  level  and  Introduced  some  uncertainty  In  the 
radio  spectra  obtained  for  these  two  coronal  holes. 

•  Chapter  5  describes  the  use  of  the  observed  radio  spectra  to  produce  models 
of  Che  solar  atmosphere  Inside  coronal  holes.  A  ray  tracing  program  Is  used  In 
conjunction  with  a  standard  model  for  Che  quiet  solar  atmosphere  to  compute 
brightness  temperatures  at  different  radio  wavelengths.  Initially  for  the  stan¬ 
dard  model  and  then  for  a  series  of  modifications  affecting  the  temperature, 
electron  density  and  neutral  particle  density  profiles  of  the  standard  model  In 
an  effort  to  find  the  profiles  that  would  reproduce  the  radio  spectra  obtained 
from  our  observations. 

The  results  of  this  analysis  Indicate  that  Inside  coronal  holes  the  follow¬ 
ing  changes  occur:  The  chromosphere  expands  in  height  by  about  a  factor  of  2, 
while  the  temperatures  and  electron  and  neutral  particle  densities  decrease  to 
about  80%  of  their  quiet  sun  values.  The  critical  transition  zone  expands  by  a 
factor  of  about  5-6,  while  the  temperatures  decrease' to  about  50%  and  the  elec¬ 
tron  densities  to  about  20%  of  their  quiet  sun  values^’  Finally,  the  corona 
shows  no  substantial  expansion,  but  the  temperatures  and  electron  densities 
decrease  to  about  80%  of  their  quiet  sun  values. 


Appendix  1  lists  the  model  for  the  quiet  solar  atmosphere  used  in  this 
study.  Appendix  2  gives  a  complete  listing  of  the  computer  program  used  In  the 
computation  of  the  radio  spectra  for  modified  atmospheric  profiles,  as  well  as 
a  sample  input  and  output. 
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INTRODUCTION 


1 .  Sc  lent  If Ic  Object  tve  s 

The  purpose  of  this  contract  was  to  study  coronal  holes 
at  radio  wavelenp, ths  and  to  use  these  radio  observations  to: 
1)  Identify  coronal  holes  In  solar  radio  maps  at  different 
wavelengths,  2)  model  the  structure  of  a  coronal  hole,  1) 
compare  different  coronal  hole-;  in  order  to  see  if  there 
were  observable  differences  at  radio  wavelengths,  and  4) 
correlate  the  coronal  hole  observations  with  high  speed 
solar  wind  streams. 

The  structure  of  coronal  holes  is  still  not  well 
understood.  Efforts  to  model  coronal  holes  have  been  made  by 
several  authors  (Wlthbroe  and  Wang,  1972;  Mariska,  1973; 
Mariska  and  Wlthbroe,  1978)  using  EUV  data.  In  addition, 
Chiuderl-Drago  (Chluderi-Drago,  et  al.,  1977;  Ch luder i-Drago 
and  Poletto,  1977;  Ch iude r i-D rago ,  1980)  has  made  several 
attempts  to  model  coronal  holes  based  on  radio  observations. 
Finally,  Dulk  and  his  co-workers  (Dulk,  et  al.,  1977)  have 
attempted  to  match  the  models  deduced  from  EUV  and  radio 
data.  Unfortunately,  the  models  derived  from  EUV  data  are  in 
substantial  disagreement  with  those  derived  from  radio  data. 
It  is  important  to  resolve  these  differences  and  much  of  our 
work  has  been  aimed  toward  this  problem. 

A  clear  understanding  of  the  transition  zone  is  of  great 
importance  to  our  understanding  of  the  coronal  hole 
phenomenon.  As  Wlthbroe  (  1978)  has  pointed  out,  "Tlie 
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transport  of  encrsy  by  thermal  conduction  Into  the 
transition  zone  Is  a  primary  coronal  cooling  mechanism." 
Thus,  the  temperature  and  particle  density  gradients  in  the 
transition  zone  are  of  critical  Importance  in  understanding 
the  energy  balance  of  a  coronal  hole. 

The  link  between  large  equatorial  coronal  holes  and  high 
speed  solar  wind  streams  observed  at  the  earth  has  been  well 
established  (Mewpert  and  V izzo ,  1976;  Nolte  et  al.,  1976, 
19  78;  ilroussard,  et  al.,  1973).  Smaller  coronal  holes  and 
coronal  holes  at  high  solar  latitudes  have  also  been  related 
to  high  speed  solar  wind  streams  (Coles  and  Pvickett,  1  97  6  ; 
Rlckett,  et  al.  1976;  Sine  and  Rickett,  1973;  Baker  and 
Papagiannls,  1980),  although  the  correlation  is  less  clear 
for  these  cases,  A  qualitative  understanding  of  the 
relationship  between  a  coronal  hole  and  high  speed  wind 
streams  emanating  from  the  coronal  hole  has  been  provided  by 
Kopp  and  Orrall  (1977),  but  a  quantitative  understanding 
requires  a  detailed  knowledge  of  the  energy  balance  in  the 
coronal  hole.  T!\us  we  are  again  brought  back  to  the 
conclusion  that  consistent  models  of  coronal  holes  are 
required  to  complete  our  understanding  of  the  solar  wind 
s  t  reams . 

This  final  report  provides  a  complete  summary  of  our 
work  on  coronal  holes,  using  radio  observations  as  well  as 
interplanetary  scintillation  and  coronal  hole  maps  made  from 
P:IIV  data,  X-Ray  data,  and  the  Helium  10830  K 
spectrohellograms  made  at  Kitt  Peak.  Chapter  1  describes  the 
Identification  of  coronal  holes  on  the  synoptic  solar  radio 
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maps  made  for  other  purposes  at  Stanford  during  the  Skylab 
period  In  1973,  as  well  as  the  determination  of  a  coronal 
hole  radio  spectrum  from  different  observations  made  during 
the  same  period.  Chapter  2  describes  the  work  we  have  done 
on  the  correlation  between  high  latitude  coronal  holes  and 
high  speed  solar  wind  streams  observed  far  above  or  below 
the  ecliptic  by  Interplanetary  scintillation  measurements. 

Chapter  3  is  devoted  to  our  observations  made  at  the 
Haystack  Radio  Observatory  and  at  the  Arecibo  Radio 
Observatory  during  1977.  The  work  described  in  Chapters  1  - 
3  has  been  described  In  detail  in  our  three  previous 
Scientific  Reports  ( AFGL-TR-77-0292  ,  AFGL-TR-73-0284 , 
AFGL-TR-30-0123) . 

Chapter  4  Is  devoted  to  the  presentation  and  analysis  of 
the  solar  radio  maps  made  during  the  month  of  July,  1978. 
These  maps  were  made  as  part  of  an  international 
co-operative  effort  organized  by  us  to  provide  solar  maps 
over  a  wide  range  of  radio  wavelengths.  Finally,  Cliapter  5 
describes  the  modeling  process  we  have  used  in  analyzing  the 
radio  spectra  we  have  obtained.  This  modeling  process 
provides  us  with  a  vertical  profile  of  temperature  and 
electron  density  inside  a  coronal  hole. 

In  addition,  there  are  two  appendices.  Appendix  1  gives 
the  electron  density,  neutral  particle  density,  and 
temperature  profiles  of  the  photosphere,  chromosphere, 
transition  zone  and  corona  that  we  have  used  as  a  standard 
model  of  the  quiet  solar  atmosphere.  Appendix  IT  gives  a 
listing  of  the  computer  program  that  was  used  to  modify  the 
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stanriaril  solar  atmosphere  model  and  then  trace  ray  paths 
through  the  atmosphere  to  determine  radio  brightness 
temperatures  at  different  wavelengths.  In  addition  to  the 
program  Itself,  we  have  Included  a  sample  Input  and  output 
of  this  p  rog  ram . 

2 .  Personnel 

The  period  covered  by  this  contract  extended  from  July 
1,  1977  to  September  30,  1980.  The  principal  Investigator 
was  !’rof.  Michael  0.  Papaglannls.  From  July,  1977  to 
September,  1978,  Dr.  Fred  Wefer  worked  on  the  project  as  a 
post-doctoral  research  associate,  while  from  January  1, 
1979  to  the  end  of  the  contract  period,  this  position  was 
held  by  Or.  Kile  Baker.  At  various  times  during  the  period 
covered  by  the  contract,  three  students  (Terry  Varner, 
Micliael  Van  Steenberg,  and  Michele  de  la  Pena)  worked  on 
various  tasks  of  the  project. 

3  .  Publ icat Ions 

During  the  period  covered  by  the  contract,  two  papers, 
four  scientific  reports  including  this  final  report,  and 
five  abstracts  have  been  published  pertaining  to  the  work  of 
the  contract.  These  publications  are  listed  below: 

Papers 

"Coronal  Holes  at  11.5  and  21  cm  observed  with  the 
Areclbo  Radio  Telescope",  M.D.  Papaglannls  and 
F.L.  u'cfer;  Mature,  273.  520  ,  1978. 


"Synoptic  Charts  of  Solar  9.1  cm  and  Coronal  Hole 
Data",  F.L.  Wefer  and  M  .  D  .  P  a  pa  I  a  n  n  i  s  ;  Solar 
Physics,  1^,  13,  1980. 

Scientific  ileports 

"The  Radio  Spectrum  of  Coronal  Hole  1",  F.L.  Uefer  and 
M . D .  Papaglannls;  AFGL-TR-  7  7-029  2  ,  Soston  Unlv., 
Ser.  II,  No.  65,  1977. 

"Synoptic  Charts  of  Coronal  Holes  and  9.1  cm  Solar  Data 
from  May  Through  August  1973",  F.L.  Wefer,  M.D. 
Papaglannls,  M . F .  Van  Steenberg,  T.M.  Varner; 
AFGL-TII-78-023A,  Boston  Unlv.,  Ser.  II,  No.  63, 
1  9  79  . 

"Correlation  of  High  Latitude  Coronal  Holes  with  Solar 
Wind  Streams  Hl^h  Above  or  Below  the  Ecliptic", 
K.B.  Baker  ana  l.D.  Papaglannls;  AFGL-TR-80-0 1 2 3 , 
Boston  Unlv.,  Ser.  II,  No.  71,  1930. 

"Radio  Studies  o  Coronal  Holes",  M.l).  Papaglannls  and 
K.B.  Baker;  A FG L-TR-3 1 -00 5 3 ,  Boston  Unlv.,  Ser. 
II,  No.  72,  1981 . 

Abs  t  rac t  s 

"Coronal  Hole  Observations  at  1420  and  2600  MHz  with 
the  Areclbo  Radio  Telescope",  M.D.  Papaglannls  and 
F.L.  Wefer;  Bull.  AAS,  9,  617,  1977. 

"Observations  of  Enhanced  Radio  Emission  at  15  C’.Hz  from 
a  Coronal  Hole  Region",  F.L.  Wefer  and  M.D. 
Papaglannls;  Dull.  AAS,  9_,  6  1  7  ,  1977  . 

"The  Radio  Spectrum  of  Coronal  Holes",  M.D. 
Papaglannnls  and  F.L.  Wefer;  EOS,  5_9^ ,  3  6  7  ,  1  973  . 


Synoptic  Charts  of  Coronal  Soles  and  9.1  cm  Data  from 
Mav  Til  r  Oil  "It  Aiiisust  1971",  K.L.  Vlcfer,  A.i, 

I’.i 'i.i",  i  .1  nil  i  s  ,  'I.h.  Van  Steeiiberi',  and  T  .  1 .  Varner; 

hull.  AAS,  j_0,  684,  1  978  . 

Solar  Wind  Streams  at  High  Solar  Latitudes  and  Their 
Correlation  to  Coronal  Moles",  M.D.  Papagiannis 
and  K.3.  Baker,  EOS,  1  095  ,  1980  . 


-  1 1  - 
Chapter  1 

KADIO  OBSERVATIONS  OF  CORONAL  HOLES  FROM  THE  SKYLAB  PERIOD 

1.1.  The  Radio  Spectrum  of  Co ronal  Hole 

Althoui;(h  by  197  7  a  number  of  groups  had  observed  coronal 
holes  at  radio  wavelengths  (Du  Ik  and  Sheridan,  1974;  Furst 
and  Hirth,  1975;  Nefer,  et  al.,  1976;  Wefer  and  Blelweiss, 
1976;  Rundu  and  Liu,  1976;  Chluderi-Drago,  1974;  Lantos 
and  Avingnon,  1975;  Covington,  1976),  no  one  had  analysed  a 
single  coronal  hole  at  many  different  wavelengths  covering  a 
wide  range  of  the  radio  spectrum.  During  the  eight  day 
period  from  June  2  5  ,  197  3  to  July  2  ,  1  973  the  coronal  liole 
designated  by  the  Skylab  Workshop  on  coronal  holes  as  CHI 
was  visible  on  the  solar  disk.  On  June  28,  1973  solar  radio 
maps  were  made  independently  by  several  groups  at  several 
different  radio  observatories  and  different  radio 
wavelengths;  3.6  inm  and  2.0  cm  -  La  Posta  Observatory,  3.8 
cm  -  NSROC  Haystack  Observatory,  9.1  cm  -  Stanford  Radio 
Astronomy  Institute,  21  cm  -  University  of  Sydney  La  Fleurs 
Radio  Astronomy  Field  Station.  On  subsequent  passages  of  the 
coronal  hole,  it  was  observed  at  three  additional 
wavelengths;  2.8  cm  -  Bonn  lOO-m  radio  telescope,  1.8  m  and 
3.75  m  -  Culgoora  radioheliograph. 

Wefer  and  Papaglannls  (1977)  analysed  these  radio  maps 
in  detail  and  presented  a  radio  spectrum  of  this  hole  by 
choosing  a  point  Inside  the  coronal  hole  and  normalizing  the 
brightness  temperatures  at  each  radio  wavelength  with  the 
corresponding  brightness  temperatures  of  the  solar  disk 
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outside  the  coronal  tiole  at  the  same  radial  distance  from 
the  center  of  the  disk.  The  combined  observations  were  used 
to  produce  the  radio  spectrum  shown  In  Figure  1.  This 
analysis  Indicated  an  enhancement  of  approximately  2%  In  the 
radio  emission  around  2.0  cm,  which  was  consistent  with  the 
previous  observations  by  Wefer,  et  al.  (1976).  At  2.8  cm  the 
bri,^;htness  temperature  of  the  coronal  hole  is  beginning  to 
be  clearly  less  tlian  the  solar  disk  value,  and  by  21  cm  the 
brightness  temperature  lias  decreased  to  about  30%  below  the 
normal  disk  value.  The  minimum  of  the  spectrum,  l.e.  the 
lowest  value  of  the  brightness  temperature  inside  the 
coronal  hole  relative  to  the  brightness  temperature  outside 
the  coronal  hole,  appears  to  occur  somewhere  between  21  cm 
and  1  m . 

The  small  enhancement  (about  2%)  at  2.0  cm  was  then  and 
still  remains  a  controversial  point.  The  enhancement  appears 
to  be  quite  ilefinlte  in  the  1973  observations.  The  question 
of  a  possible  enhancement  at  2.0  cm  in  other  coronal  holes 
will  be  discussed  in  more  detail  in  Chapter  4  when  we 
consider  the  1978  observations. 

One  other  imoortant  feature  of  our  analysis  of  CUl  radio 
data  was  tlie  discovery  that  the  coronal  hole  could  indeed  be 
observed  in  the  9.1  cm  radio  maps  from  Stanford.  The 
Stanford  people  themselves  had  claimed  to  be  unable  to  see 
coronal  holes  in  their  data.  The  fact  that  our  radio 
spectrum  indicates  that  one  ought  to  be  able  to  see  coronal 
holes  around  9  cm  as  10-20%  depressions  in  the  brightness 
temperature  relative  to  the  solar  disk  background  lead  us  to 
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i;o  back,  tbroui’h  all  the  9.1  cm  Stanford  maps  that  were 
available  durini;  the  Skylab  period  in  a  search  for  coronal 
t\oles.  The  results  of  tliat  search  are  presented  in  the  next 
sect  ion . 

A  detailed  description  of  our  work  on  CHl  is  available 
in  our  Scientific  Report  Mo.  1,  AFGL-TR- 7 7-029  2  ,  Oec.  1977 
by  K.L.  Wefer  and  T.D.  Papagiannis. 


1.2.  Identification  o  f  Coronal  il  o  1  e  s  on  the  Stanford  9 . 1  cm 
Synoptic  Charts 

The  period  from  !tay  28  ,  1973  to  August  1  1  ,  1  973  was  an 
ideal  period  for  the  study  of  coronal  holes  because  it  was 
the  only  time  at  which  solar  X-Ray  data  from  ATM-Skylab  and 
9.1  cm  radio  maps  from  the  Stanford  Radio  Astronomy 
Institute  wore  simultaneously  available.  Tn  addition,  it  was 
a  period  near  solar  minimum  when  large  transequatorial 
coron.il  iiolcs  were  present  on  the  sun. 

Our  study  of  the  radio  spectrum  of  coronal  hole  Cill 
('.v'efer  and  Pap  a  g  Iannis,  1977)  had  demonstrated  that  a 
coronal  hole  could  easily  be  seen  at  9.1  cm.  Graf  and 
liracewell  (1975)  liad  produced  synoptic  charts  of  the  9.1  cm 
data  for  each  solar  rotation  from  June,  1962  up  to  August, 
1973,  when  the  9.1  cm  system  was  shut  down.  In  their  charts 
the  emphasis  was  placed  on  active  regions  and  the  lowest 
contoured  level  was  at  a  brightness  temperature  of  40,000  K. 
The  central  disk  brightness  temperature  for  the  quiet  sun  at 


9.1  cm  Is  approximately  30,000  K,  and  since  at  this 
wavelength  coronal  holes  manifest  themselves  as  10-20% 
depressions  in  the  brightness  temperature  relative  to  the 
quiet  sun  background,  it  is  understandable  why  the  synoptic 
charts  produced  by  Oraf  and  'iracewell  did  not  show  any 
coronal  holes. 

Graf  and  Bracewell  kindly  provided  us  with  the  computer 
tapes  of  the  dally  data  of  their  9.1  cm  observations,  and 
from  these  data,  we  produced  our  own  synoptic  charts, 
emphasizing  contours  near  and  below  the  quiet  sun  level.  The 
basic  procedure  was  a  simple  one.  On  each  day,  we  extracted 
the  data  corresponding  to  the  central  '.iieridlan  and  collected 
30  days  worth  of  data  to  form  a  grid  from  which  a  synoptic 
contour  map  was  drawn  covering  slightly  more  than  one 
Carrington  rotation.  Since  this  method  produced  charts  based 
only  on  central  meridian  data,  they  are  suitable  only  for 
studying  long-lived  coronal  holes.  As  a  check  on  the  resvilts 
obtained  in  this  fashion,  we  also  constructed  a  few  charts 
using  data  just  to  the  left  or  right  of  t)ie  central 
meridian.  These  charts  showed  the  same  general  features  as 
the  central  meridian  charts,  but  differed  somewhat  in  the 
finer  (and  shorter-lived)  details. 

A  typical  synoptic  chart  Is  shown  In  Figure  2.  The 
shaded  areas  show  the  positions  of  coronal  holes  as 
determined  from  the  ATM-Skylab  data  (Molte,  et  al.,  1976b) 
for  rotation  1602.  In  this  chart,  CMl  (near  longitude  zero) 
and  Ctl3  (near  longitude  260°)  are  clearly  manifested  as 
depressions  in  the  9.1  cm  emission.  Coronal  hole  013,  as 
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well  as  several  smaller  coronal  holes,  did  not  show  up  in 
the  9.1  cm  chart,  probably  because  their  sizes  are  smaller 
than  the  beamwidth  of  the  Stanford  antenna  and  the 
additional  side-lobe  interference  from  nearby  active 
r eg i ons  . 

A  complete  description  of  this  work  is  Included  in  our 
Scientific  Report  No.  2,  AFCL-TR-78-0284 ,  January  1979  by 
F.L.  Wefer,  M.D.  Papagiannis,  T.  Varner  and  ?i.  Van 
Steenberg,  and  in  the  paper,  “Synoptic  Charts  of  Solar  9.1 
cm  and  Coronal  Hole  Data"  by  F.L.  Wefer  and  M.D.  Papagiannis 
(1980). 
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Chapter  2 

HIGH  CATITUUK  CORONAL  HOLRS  ANO  TMEIR  CORRELATION  TO 
MTG!I  LATITUDE  INTERPLANETARY  SCINTILLATION 

2.1.  Xlie  Correlation  o  f  IPS  and  Coronal  Mole  Data 

Althoiip.h  the  correlation  bctia^cen  large  equatorial 
coronal  holes  and  high  speed  solar  wind  streams  observed 
near  the  earth  has  been  well  established  (Krieger,  et  al., 
1971;  Kopp  and  Orrall,  1977;  Mundhausen,  1977),  the 
generalization  that  all  coronal  holes  are  sources  of  high 
speed  solar  wind  streams  is  not  well  established.  In  fact, 
Mundhausen  (1977)  has  suggested  that  the  correlation  may  not 
hold  true  for  smaller  coronal  holes. 

A  general  correlation  between  large  coronal  holes  and 
high  solar  wind  velocities,  derived  from  interplanetary 
scintillation  (IPS)  measurements,  was  obtained  in  the  70's 
at  the  University  of  California,  San  Diego  by  Rickett,  et 
al.  (1976).  They  found  that  about  half  of  the  large  coronal 
hole  passages  seen  over  several  solar  rotations  were 
associated  with  high  speed  solar  wind  streams. 

Recently,  Coles,  et  al.  (1980)  have  correlated  the 
average  solar  wind  speed  observed  at  high  solar  latitudes 
with  the  size  of  the  polar  coronal  hole.  They  have  shown 
that  as  the  solar  cycle  approaches  maximum,  the  polar 
coronal  hole  decreases  in  size  and  the  average  solar  wind 
speed  at  high  solar  latitudes  diminishes. 

Nevertheless,  a  definitive  correlation  showing  a 
connection  between  specific  high  latitude  coronal  holes  and 


high  speed  solar  wind  streams  corresponding  to  similar  hl/.h 

solar  latitudes  had  not  been  attempted.  To  this  end  we  made 

a  more  detailed  correlation  analysis  between  the  velocities 

of  solar  wind  streams  high  above  or  below  the  ecliptic 

obtained  from  IPS  data,  and  the  positions  and  shapes  of  lilgh 

solar  latitude  coronal  holes  obtained  from  coronal  hole 

o 

synoptic  maps  constructed  from  the  1  0  3  30  A  Hellu-i 

spectrohellograms. 

For  the  past  several  years  there  have  been  only  two 
groups  making  dally  interplanetary  scintillation 
measurements,  namely,  the  group  headed  by  Coles  at  UCSD  and 
the  group  headed  by  Watanabe  and  Kaklnuiaa  at  the  Toyokawa 
Observatory  In  Japan.  We  utilized  both  sets  of  data  In  our 
correlation  In  order  to  provide  as  complete  a  coverage  in 
time  as  possible  and  whenever  possible  to  provide  a  clieck  of 
one  set  against  the  other.  Our  study  covered  the  two  and  a 
half  year  period  from  January  I,  1977  to  June  30,1979.  For 
the  first  year  and  a  half  (Jan. I,  1977  to  June  30,1973)  we 
had  IPS  data  from  both  UCSI)  and  Toyokawa.  For  the  last  year, 
however,  we  only  had  data  from  UCSD  because  the  Toyokawa 
data  had  not  yet  been  processed  and  published. 

The  position  and  shapes  of  coronal  holes  were  determined 
from  the  published  monthly  synoptic  charts  (Solar 
Geophysical  Data  Prompt  Keports)  made  by  Dr.  J.  Harvey  from 
from  the  Kltt  Peak  Helium  10830  %  spectrohellograms.  It 
should  be  noted,  however,  that  because  of  the  difficulties 
Involved  In  interpreting  these  spectrohellograms,  there  is 
some  uncertainty  in  the  details  and  exact  boundaries  of  the 


coronal  holes  shown  In  the  monthly  synoptic  charts.  In 
addition,  for  the  coronal  holes  present  in  July,  1978,  wo 
have  also  made  use  of  our  radio  observations  to  determine 
the  structure  and  boundaries  of  the  coronal  holes.  This  work 
is  discussed  in  detail  in  Chapter  U  and  the  radio  maps  are 
shown  in  KlRures  23  -  44. 

The  method  of  analysis  was  a  simple  one.  If  the  point  of 
maximum  scintillation,  P,  for  a  particular  radio  source  (see 
Figure  3)  corresponded  to  a  solar  latitude  that  would 
Intersect  the  boundaries  of  a  coronal  hole,  we  looked  at  the 
IPS  data  to  determine  if  a  high  speed  solar  wind  stream  was 
present.  We  would  expect  the  coronal  hole  to  produce  a 
scintillation  effect  (a  high  speed  wind  stream)  when  it  is 
pointed  toward  tlie  point  of  maximum  scintillation.  The 
actual  effect,  however,  would  be  delayed  by  approximately 
one  day,  which  is  tlie  time  it  would  take  for  the  hiv'd  soeeJ 
./ind  stroan  to  advance  from  tlie  coronal  hole  to  the 
scintillation  point.  The  actual  magnitude  of  the  delay  Is  a 
function  of  the  solar  wind  speed  and  the  distance  of  the 
scintillation  point  from  the  sun  and  can  be  computed  through 
an  iterative  procedure. 

To  determine  the  correlation  we  had  to  define  a  set  of 
criteria  for  determining  how  good  a  particular  correlation 
was.  We  defined  four  categories  for  the  correlation; 

1)  Very  Good  -  High  speed  solar  wind  streams  were 
observed  on  the  expected  dates.  The  increase  in 
wind  speed  was  significant  (25%  or  more).  If  Che 
boundary  of  the  coronal  hole  was  such  that  Che 
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scintillation  point  passed  in  and  out  of  the 
coronal  hole  area  more  than  once,  this  should  be 
reflected  in  the  IPS  data. 

2)  Good  -  A  significant  increase  in  the  solar  wind, 
followed  by  a  corresponding  decrease,  must  have 
occurred  within  2  days  of  the  expected  dates. 

3)  Fair  -  The  solar  wind  speed  must  liave  been  higher 
than  average  during  most  of  the  expected  period 
(with  a  tolerance  of  2  days). 

4)  Poor  -  The  solar  wind  speed  stayed  near  the 
average  value  or  decreased  during  the  dates  of 
expected  high  spjeed  solar  wind. 

Utilizing  t\\ese  criteria,  we  were  able  to  analyse  a  total  of 
26  cases  and  classify  them  in  one  of  the  above  four 
categories.  We  were  unable  to  classify  seven  other  cases  due 
to  a  lack  of  all  the  needed  data.  As  can  be  seen  from  the 
hlstog.ram  of  the  26  cases  analysed  (Figure  4),  the 
correlation  in  about  ^07,  of  the  cases  was  either  "Very  Good" 
or  "Good",  while  in  fewer  than  25X  of  the  cases  the 
correlation  was  "Poor". 

Considering  the  assumptions  (radial  propagation  of 
solar  wind  streams,  constant  solar  wind  velocity  along  tlie 
radial  path,  etc.)  Involved  in  this  analysis,  the  sizable 
uncertainties  In  the  solar  wind  speeds  deduced  from  IPS 
i.ieasuroments ,  as  well  as  the  uncertainties  in  the  positions 
and  boundaries  of  coronal  holes  deduced  from  the  10830  A 
synoptic  charts,  the  overall  correlation  is  surprisingly 
.’ood.  We  conclude  that  It  is  certainly  true  that  a 


solar 
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s  i  ■ n  1  f  i  c  a  n  t  fraction  (probably  most)  of  the  high 
latitude  coronal  holes  are  sources  of  high  speed  solar  wind 
streams.  N'evertheless,  It  is  also  quite  possible  that,  like 
the  equatorial  coronal  holes,  some  high  latitude  coronal 
holes  do  not  produce  high  speed  solar  wind  streams  on  a 
continuous  basis,  although  at  various  stages  in  their 
evolution,  they  may  do  so.  It  is  also  quite  evident  that  the 
strength  of  the  solar  wind  streams  is  related  to  the  size  of 
the  coronal  hole.  It  is  quite  possible  therefore  tliat  some 
of  the  smaller  streams  are  totally  missed,  either  because  of 
tbelr  relatively  sliort  duration  and  our  limited  means  of 
observation,  or  because  they  get  dissipated  or  diverted  in 
interplanetary  space  before  they  reach  the  expected  point  P 
of  maximum  scintillation. 

A  complete  description  of  tliis  study,  along  with  the 
data  used  is  given  in  our  Scientific  Report  No.  3, 
AKGL-TR--i()-0123,  April,  1980  by  K .  8 .  Raker  and  M.D. 
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Chapter  3 

OBSERVATIONS  OF  CORONAL  HOLES  AT  ARECIBO  AND  HAYSTACK 

IN  SEPTEMBER  OF  1977 

3.1.  Coronal  Hole  Obse  rva  t  tons  a  t  11.5  and  2 1  cm  with  the 
A r e c 1 b o  Radio  Telescope  1 n  19  7  7 

During  the  five  day  period  from  August  31  to  September 
4,  1977,  we  obtained  seven  maps  of  the  sun  at  11.5  cm  and 
21  cm  using  the  Arecibo  radio  telescope.  These  solar  maps 
were  made  by  two  different  methods.  In  the  drift  scan  method 
the  telescope  remained  stationary  and  the  sun  was  allowed  to 
drift  past  the  beam  of  the  antenna.  Each  scan  took  about  4 
minutes  and  the  total  map,  consisting  of  21  scans,  took 
about  1.5  hours.  In  the  second  method  the  antenna  was  driven 
at  twice  the  siderial  rate  while  executing  a  boust rophedonlc 
raster.  This  method  made  it  possible  to  cut  the  time 
required  for  each  map  in  half  and  allowed  us  to  make  naps  at 
both  11.5  cm  and  21  cm  during  the  approximately  two  hour 
period  around  local  noon  that  the  sun  can  be  observed  with 
the  Arecibo  radio  telescope  at  this  time  of  the  year.  We 
were  thus  able  to  make  seven  solar  maps  during  this  five  day 
observing  period. 

This  was  the  first  time  that  the  Arecibo  radio 
telescope  was  used  to  make  solar  maps.  The  Arecibo  dish  was 
of  unique  value  for  this  project  because  of  Its  high 
sensitivity  (“1%  of  solar  disk  background)  and  its  high 
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aiu;ul.ir  resolution  (1-4  arcmln)  at  decimeter  (11.5  and  21 
cn)  wavelen;;  tlis.  The  discovery  that  coronal  lioles  become 
clearly  visible  at  these  wavelengths  confirmed  the 
lin[iortance  of  the  Arecibo  radio  telescope  for  decimetrlc 
solar  observations. 

The  results  of  this  project  are  shown  in  Figure  5,  which 

shows  five  of  the  solar  maps  made  during  this  five  day 

period.  The  coronal  hole  boundaries  Indicated  on  each  map 

o 

were  obtained  from  the  Kltt  Peak  10830  A  spectroheliograms, 
which  wore  provided  by  Ur.  J. Harvey  of  the  Kitt  Peak 
Observatory.  Tliese  spectroheliograms  are  used  for  an  optical 
identification  of  coronal  holes  in  the  absence  of  X-ray  or 

O 

UV  observations.  The  10830  A  spectroheliogram  for  September 
2,  l'^77  is  also  siiown  in  Figure  5. 

Tlie  coronal  hole  in  the  southern  hemisphere  is  very 
obvious  and  is  very  similar  in  shape  and  position  to  the 

O 

coronal  hole  seen  in  the  10830  A  s  pe  c  t  roll  e  1 1  og  ram  .  The 
coronal  liole  in  the  northern  liemlsphere  was  much  narrower 
anti  was  therefore  more  difficult  to  observe  at  radio 
wa  ve  1  •;  ng  t  h  s  .  Nevertheless,  a  depression  in  radio  emission  is 
clearly  visible  in  the  region  of  the  north  pole. 
L'nfortiinately,  the  picture  in  the  north  is  confused  by  the 
presence  of  several  filaments  (indicated  by  dark  wiggly 
lines)  which  also  show  up  as  depressions  at  radio 
wavelengths  (Straka,  Papagiannis  and  Kogut,  l')75). 

This  work  Is  more  fully  reported  in  the  article, 
"Coronal  holes  at  11.5  and  21  cn  observed  with  the  Arecibo 
radio  telescope",  in  Nature  by  M.U.  Papagiannis  and  F  L. 
W  e  f  e  r ,  19  7  8  . 


.  2  .  C o  r  ()  n a  I  h  e  r  v  a  i  ions  ■>  t  II a  y  s  t  .i  c  k  a  nil  \.n  P o  s  t  a  i  n 
197  7 

In  addition  to  the  Arectbo  radio  naps,  we  were  also  able 
diirini^  the  same  period  to  make  radio  maps  at  3.81  cm  and 
1.9  4  cm  usint;  the  Haystack  Observatory  radio  telescope. 
These  maps,  alon';  with  maps  at  2.0  cm  and  0.86  cm  from  the 
La  Posta  Observatory  are  shown  in  Figure  6.  The  position  of 
the  coronal  holes  and  filaments  are  also  indicated  on  these 
maps.  The  coronal  holes  appear  as  slight  depressions  at  3.31 
cm.  At  2.0  cm  and  1.94  cm,  the  holes  appear  as  slight 
enhancements.  This  apparent  enhancement  will  be  discussed  in 
more  detail  in  Chapters  4  and  6. 

3.3.  The  Rad lo  Spec  t  rum  o  f  the  September  1  977  Coronal  Hole 

From  our  observations  at  Areclbo  and  Haystack,  combined 
with  additional  solar  radio  mans  made  by  tlie  La  Posta 
observatory  in  California,  we  were  able  to  determine  a 
partial  radio  spectrum  of  this  coronal  hole  (Figure  7).  The 
results  are  quite  consistent  v^ith  those  for  CHI  in  1973 
(Figure  1  in  Chapter  1)  particularly  with  respect  to  the 
approximately  2%  enhancement  around  2  cm. 
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FIGURE  7 


Chapter  4 


CORONAL  HOLK  RADIO  OBSERVATIONS  DURING  JULY  197U 

4,1.  International  Effort 

From  the  previous  work  it  had  become  clear  that 
simultaneous  observations  of  coronal  holes  over  a  wide  r.in;>,e 
of  radio  wavelenp,  ths  would  allow  us  to  determine  the  radio 
spectr<im  of  coronal  holes,  which  would  be  of  ^rcat  value  in 
our  effort  to  understand  the  Internal  structure  and  the 
pliyslcs  of  coronal  holes.  It  would  also  he  useful  to  have  a 
certain  amount  of  redundancy  in  the  observations  to  insure  a 
higher  level  of  dependability  on  the  data. 

Starting  early  in  1978  we  began  to  organize  an 
international  co-operative  effort  to  observe  coronal  holes 
at  many  different  radio  wavelengths.  The  project  was  well 
received  by  many  groups  around  the  world  and  the  tine  period 
selected  for  this  combined  effort  'vas  July  of  197  3.  The 
month  of  July  was  chosen  primarily  because  we  wanted  to 
include  observations  at  21  cm  and  11.5  cm  from  the  Arecibo 
observatory  which  can  observe  the  sun  only  during  the  summer 
and  allows  the  longest  observing  time  near  the  summer 
solstice.  In  addition,  most  of  the  participating 
observatories  were  located  in  the  northern  hemisphere  and 
would  be  expected  to  have  better  observing  conditions  during 
the  summer  months.  In  order  to  Increase  the  probability  that 
a  coronal  hole  would  be  present  on  the  solar  disk  Jurinv,  our 
observations,  wc  decided  to  use  two  separate  observing 
periods  of  about  five  days  each  separated  by  about  ten  days. 
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The  two  periods  were  July 

Observatory  Location 

6-10  ; 

TABLK 

Size 

and  July  15-1 3 , 

1 
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s  iiinain  r  i  ze  fl  la  Table  1.  Unfortunately  the  radio  maps  from  the 
Tovokawa  observatory  and  the  University  of  British  Columbia 
had  severe  slJelobe  problems  and  difficulties  with  the 
processing  of  their  radio  data  which  made  these  maps 
virtually  useless  for  studying  coronal  holes.  These  maps, 
therefore,  were  not  Included  in  our  analysis. 

4.2.  First  Period  Obse  rva  1 1  ons  (July  6-  10  ,  1978) 

The  first  period  chosen  for  tlie  radio  observations  was 
from  July  6  to  July  10,  1978.  Not  all  observatories  were 
able  to  make  observations  for  the  entire  period,  but  there 
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Is  nonetheless  ;i  reasonable  overlap  of  dates  from  all 
observatories.  Table  2  1 Ists  the  solar  parameters,  P,  Bo  , 
and  Lq  (Carrington  Longitude  of  the  Central  ‘lerldian)  at 
0000  UT  for  each  of  the  six  days  of  the  first  period. 


TABLE  2 

Solar  Parameters  at  0000  UT 


DATE 

la 

P 

7-6 

23  5 

+  3.4 

-0.6 

7-7 

272 

3 . 5 

-0.1 

7-3 

259 

3.6 

+  0.3 

7-9 

246 

3.7 

+0.3 

7-10 

232 

3.8 

+  1  .2 

Table  3  lists  the  participation  observatories  and  indicates 
the  wavelengths  and  the  days  on  which  each  observatory  made 
its  observa  t ions . 


TABLE  3 

Participating  Observatories  -  First  Period  Observations 


July  Date 


Observa  tory 

Wavelength(cm)  6 

7 

8 

9 

10 

La  Pos  t  a 

0.86 

X 

X 

X 

X 

X 

La  P 0 s  t  a 

2.00 

X 

X 

X 

X 

'lay  stack 

0.69 

X 

X 

day  stack 

1  .38 

X 

X 

X 

Kay  stack 

1.93 

X 

X 

X 

Haystack 

3.31 

X 

X 

itapet inga 

0.7  0 

X 

X 

X 

I tapet Inga 

1  .39 

X 

X 

X 

A  r  e  c  1  b  o 

11.5 

X 

X 

X 

A  r e  c i bo 

2  1  .4 

X 

X 

C  u  1  g  o  o  r  a 

187 

X 

Culgoora 

375 

X 

X 

As  is  clear 

f  rom  Table  3 

,  July  8 

in  the 

date 

on 

which 

most  complete 

coverage  is 

available 

• 

The  first  period  turned  out  to  be  a  poor  one  for  the 
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purpose  of  observint;  coronal  holes.  The  Helium  10!330  A  Jata 
show  only  the  usual  small  polar  depressions.  The  helium 
synoptic  map  from  Kitt  Peak  for  half  a  Carrin;;ton  rotation 
centered  approximately  on  July  7  is  shown  in  Kigure  3.  Tliere 
are  no  coronal  holes  hut  there  are  several  active  regions 
present.  These  active  regions  can  also  be  clearly  seen  in 
the  Ha  and  Calcium  pictures  of  tlie  sun  for  tills  period.  The 
ila  picture  from  Sacramento  Peak  and  the  MacMath-Hulbert 
Calcium  Report  for  July  8  are  also  shown  in  Figure  3.  The 
three  largest  active  regions  are  McMath  15335,  15391,  and 
15393.  The  Ha  picture  shows  a  very  clearly  defined  filament 
in  the  northern  hemisphere  at  an  approximate  Carrington 
longitude  of  240P  .  Thus,  during  the  first  period  this 
filament  is  always  on  the  east  side  of  the  central  meridian. 
It  is  important  to  note  the  position  of  this  filament, 
since,  as  mentioned  above,  large  filaments  appear  very 
similar  to  coronal  holes  on  radio  maps. 


Radio  Maps  for  July  6-10, 1978 

Figures  9  -  20  show  representative  maps  from  the  first 
period  observations.  Although  some  details  vary  from  day  to 
day,  the  important  features  are  relatively  constant  over  the 
five  day  period  from  July  6  to  July  10. 

The  shortest  wavelength  data  available  are  the  0.69  cm 
maps  from  the  Haystack  observatory.  At  this  very  short 
wavelength  all  the  information  comes  from  deep  in  ttie 
chromosphere  and  all  solar  features,  including  the  active 
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rt'",  i  ons  ,  sho'.j  very  little  contrast  with  Che  siir  roiiiul  i  n;; 
Tuict  area.  Still,  tlic  active  re  p.  Ions  can  be  seen,  as  well 
as  the  su,m;estion  of  a  sli(>ht  depression  corresponding  to 
the  position  of  the  filament  in  the  northeast.  The 
It  ape  ting a  maps  at  0.7  0  cm  for  the  first  period  were  marred 
bv  spurious  readings.  We  have  cleaned  up  the  map  by 
interp(-)lation  to  remove  the  spurious  data,  but  the  resultant 
maps  show  few  features. 

The  'laystick  I.IS  cm  maps  and  the  Itapetinga  1.39  cm 
naps  noth  clearly  show  the  active  regions  and  the  filament 
in  tiie  northeast.  The  Itapetinga  maps  were  again  marred  by 
spurious  data,  but  in  regions  where  we  can  be  sure  the 
Itapetinga  naps  are  reasonably  clean,  the  features  are  quite 
consistent  with  the  Haystack  maps. 

The  La  Posta  2.0  cm  maps  and  the  Haystack  1.93  cm  maps 
are  ilso  quite  consistent.  Hoth  show  the  active  regions  as 
stroi.;  enhancements  and  the  filament  as  a  slight  depression, 
rill'  laystack  maps  also  have  a  suggestion  of  a  depression 
slightly  east  of  the  central  meridian  in  the  far  south.  The 
Oac  I'eaK  'la  map  indicates  the  presence  of  a  small  filament 
In  this  region,  which  Is  probably  the  source  of  this 
depression. 

The  Arecibo  maps  at  11.5  cm  and  21.4  cm  show  the  active 
rej;lon.s  and  the  filament  In  the  northeast  quite  clearly. 

The  long  wavelength  maps  from  the  Culgoora  Observatory 
(  1  . 'i  7  m  and  3.75  m)  show  much  less  detail  than  the  other 
radio  maps.  The  active  regions  can  be  seen  at  both 
wavelengths,  but  at  tt\e  shorter  wavelength  the  maps  are 
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confused  by  negative  sidelobes  from  these  active  regions.  \n 
additional  complication  in  interpreting  these  maps  is  that 
tlie  sun  was  very  low  in  the  sky  during  those  J<ily 
observations  at  Culgoora.  Refraction  of  these  long 
wavelength  radio  waves  through  the  earth's  atmosphere 
displaces  the  radio  image  of  the  sun  from  the  ephemerJs 
position,  making  it  difficult  to  determine  tlie  precise 
position  of  the  observed  features  on  the  solar  disk. 
Unfortunately  these  maps  were  only  available  in  the  form  of 
finished  contour  maps  and  not  In  the  form  of  data  points. 
This  made  It  impossible  for  us  to  re-contour  the  maps  w  r  t  li 
enpliasls  on  contours  near  the  quiet  sun  level. 

4.3.  Second  Period  Obse  rva  t  tons  (July  15-18,  1978  ) 

The  coverage  of  the  radio  spectrum  during  the  second 
period  was  more  complete  than  for  the  first  period.  Table  4 
lists  the  solar  parameters  for  the  four  days  of  the  second 
period  and  Table  5  lists  the  participating  observatories, 
wavelengths,  and  the  days  each  observatory  made  its 
observations.  July  16  and  17  are  the  dates  with  the  best 
c  overage . 

TABLE  4 

Solar  Parameters  at  0000  UT 


DATE 

ko 

“o 

P 

7-15 

166 

4.3 

3.4 

7-16 

1  53 

4.4 

3.9 

7-17 

140 

4.5 

4.4 

7-18 

127 

4.6 

4.3 

The  helium  10330  A  spect  rohel  logram  for  July  16  Is  sliown 
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in  Ki»,iirc  2  1,  .ind  the  Helium  synoptic  chnrt  for  the  hnlf 
0,irri.i‘>,  tion  rotation  centered  approximately  on  July  16  is 
shown  in  Kif’ure  22.  The  outlines  of  the  coronal  holes  are 
also  shown  as  they  were  drawn  by  Jack  Harvey,  of  the  Kitt 
Peak  Observatory,  who  kindly  provided  us  with  these  maps. 
Tlie  Helium  data  show  a  larp,e  coronal  hole  extending  from  the 
south  pole  towac<l  the  equator,  a  small  hole  at  latitude  1  0°  S 
and  a  somewhat  larger  hole  at  latitude  40° K.  \  number  of 
active  r  e  '>  1  o  n  s  are  also  present. 


TA3Lr,  5 


Participatin',  Observatories  -  Second  Period  Observations 


Observatory 

viavc  length(cm) 

15 

July 

16 

Date 

17 

13 

1.  a  Post  a 

0.8  6 

X 

X 

X 

X 

I.  a  Post  a 

2.00 

X 

X 

X 

X 

1  a  y  s  t  d  c  k 

0.69 

X 

X 

X 

'  1  a  y  s  t  a  <•  k 

1  .  )3 

X 

X 

X 

1  a  y  s  t  a  c  k 

1.93 

X 

X 

X 

d  a  V  s  t  a  c  k 

3.8  1 

X 

X 

X 

t  t  a  p  0  t  i  n  g  a 

0.70 

X 

X 

X 

T  t  a  n  e  t  i  n  '  a 

1  .39 

X 

X 

X 

A  r  e  c  i  h  o 

11.5 

X 

X 

X 

A  r  e  c  L  1)0 

Cu  1  >■  o  o  r  a 

2  1.4 

137 

X 

X 

X 

Y 

X 

(’.  u  1  '  o  ora 

37  5 

X 

X 

X 

The  Sac  Peak  da 

picture 

(  F  i  H  u  r  e 

22) 

for 

July 

16  also 

'lows  several  strong 

act  i  V  e 

regions 

a  s 

well 

as 

several 

filaments.  Also  sliown  in  Figure  22  is  the  McMa  t  h-!iu  1  be  r  t 
Calcium  Heport.  The  most  prominent  active  region  is  tlc'tath 
I  540  1 . 

Tl»e  coronal  hole  in  the  south  lies  in  a  relatively 
quiet  area  of  the  sun  which  should  make  it  easily  detectable 


-Sl¬ 
at  radio  wavolcn«tli.s.  On  the  other  hand,  the  coronal  hole  In 
the  northern  heinisphcre  lies  very  close  to  the  large 
filament  between  Mc'lath  15410  and  McMath  15415,  which  makes 
It  difficult  to  Identify  the  coronal  hole  at  radio 
wave le ng t  hs  . 


I 


i 


i 
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It i  o  .'In  p  s  for  .July  1  'i  -  1  8  ,  1  *)  7  8 

Kl;;iires  28  -  44  show  some  of  the  maps  for  the  second 
period  observations.  Special  emphasis  has  been  placed  on  the 
maps  for  the  two  days  July  16  and  July  17,  since  these  two 
Jays  represent  the  days  for  which  wo  have  the  most  complete 
wavelength  coverage. 

As  w  1  t  li  the  first  period  maps,  the  shortest  wavelength 
maps  show  little  contrast.  The  Haystack  0.69  era  maps  are 
virtually  devoid  of  features  other  than  the  strong  active 
regions.  Tlie  Itapetlnga  maps  at  0.70  cm  show  a  distinct 
depression  in  the  southern  hemisphere  corresponding  in 
position  to  the  coronal  liole.  In  contrast,  however,  the  0.36 
cm  maps  from  La  I’osta  show  only  tl»e  filament  and  coronal 
fiolc  in  the  nortlieast  a.s  depressions. 

Tlie  Haystack  maps  at  1.33  cm  also  show  distinct 
differences  from  the  Itapetinga  maps  at  1.39  cm.  The 
Haystack  maps  show  clear  depressions  in  both  the  north  and 
the  south.  In  both  cases,  however,  the  position  of  the 
depression  is  shifted  to  the  west  of  the  position  of  the 
coronal  hole  as  determined  from  the  Helium  10330  X  charts. 
This  apparent  shift  In  the  position  had  also  been  noted  in 
previous  observations  (We  for  and  Papaglannls,  1977). 

The  Haystack  maps  also  show  a  depression  on  the  east  limb 


in  the 

southern  hemisphere. 

The 

Ha 

pictures  show 

seve  ra 1 

small 

filaments  in  this  regi 

on  of 

the  east  limb  which 

a  re 

probably  unresolved  In  tlie 

rad  i 

O 

maps  and  sitow 

up  as 

a 

single 

extended  depression. 

The 

Itapetlnga  maps  at 

1  .39 

c  m 

show 

only  a  sliglit  suggestion 

of 

the  depression 

i  n 

t  he 
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north,  while  the  depression  in  the  south  appears  to  be 
merely  an  extension  of  the  large  depression  coming  from  the 
east  limb.  Presumably  the  differences  between  the  Haystack 
and  Itapetlnga  maps  reflects  the  difference  in  angular 
resolution  of  the  two  antennas.  At  Haystack  the  resolution 
is  about  1.6  aremin  at  1.38  cm  while  at  Itapetlnga  the 
resolution  is  only  about  4  aremin.  Thus,  the  nearby  active 
regions  mask  the  coronal  holes  in  the  Itapetlnga  maps, 
whereas  the  Haystack  antenna  was  able  to  resolve  these  as 
separate  features. 

Tlie  :iiaps  at  1.9  3  c  i  (Haystack)  and  2.0  c.i  (La  i’osta) 
noth  show  slight  depressions  corresponding  to  the  coronal 
hole  in  the  south.  As  with  the  1.38  cm  maps,  the  position  of 
the  coronal  hole  appears  to  be  shifted  to  the  west  of  its 
position  in  the  Helium  10830  X  maps.  Indeed,  this  westward 
shift  is  a  consistent  feature  of  the  longer  wavelength  maps. 
In  the  north,  the  Haystack  maps  show  the  coronal  hole  as  a 
distinct  depression.  Slightly  west  of  the  depression, 
however,  is  a  small  region  of  enhanced  emission.  This 
enhancement  does  not  correspond  to  any  clear  features  in  the 
Ha  or  Calcium  s pe c t r ohe 1 1 og rams ,  but  it  does  appear  to  be 
r  eal . 

The  La  Posta  2.0  cm  maps  Indicate  a  slight  enhancement 
in  the  region  of  the  northern  coronal  hole.  This 
disagreement  between  Haystack  maps  and  La  Posta  maps  may  be 
the  result  of  differences  In  angular  resolution.  The  regions 
of  enhanced  emission  seen  in  the  Haystack  maps  as  separate 
features  are  probably  masking  the  slight  depression  of  the 


coronal  hole  in  the  La  I'osta  maps.  If  slight  enhancements 
near  the  boundary  of  a  coronal  hole  are  common  features, 
this  may  explain  the  previous  results  of  Wefer,  et  al. 
(197b),  Wefer  and  !lleiweiss  (1976),  and  our  own  results  from 
1977  (see  Chapter  3),  where  coronal  holes  appeared  as 
e  nil  a  nc  e  me  n  t  s  of  about  I  -  2  ^  at  2.0  cm. 

The  Haystack  maps  at  3.31  cm  show  the  usual  active 
regions  as  well  as  the  coronal  holes  in  both  the  north  and 
the  south.  As  with  the  1.93  cm  maps,  these  maps  show  a  small 
region  of  enhanced  emission  slightly  to  the  west  of  the 
northern  coronal  hole. 

The  Arecibo  maps  at  both  21.4  cm  and  11.5  cm  show  both 
coronal  holes  quite  clearly,  including  the  eastward 
extension  of  the  southern  coronal  hole.  The  depression  on 
the  east  limb  is  also  clearly  visible  as  a  separate  feature. 
As  before,  the  active  regions  are  also  prominent  features. 

As  was  tlie  case  during  the  first  period  observations, 
the  Cuigoora  maps  at  1.87  m  and  3.75  m  are  more  difficult  to 
interpret  than  the  other  maps,  due  to  the  uncertainty  about 
the  position  of  the  solar  disk  and  the  coarse  contour 
intervals  in  the  solar  maps  provided.  At  3.75  m,  the  most 
prominent  feature  is  the  large  brig  lit  region  near  the 
center,  which  Is  probably  the  strong  active  region  McMath 
154  0  3.  South  of  tills  feature  is  a  slight  depression  which 
corresponds  in  position  to  the  coronal  hole.  As  usual,  the 
1.87  m  maps  show  many  more  features  than  the  3.75  m  maps. 
Again,  sldelobcs  create  major  difficulties  in  interpreting 
these  naps.  There  appear  to  be  two  regions  of  depressed 
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emission  In  the  south,  with  the  eastern  one  corresponding 
more  closely  In  position  to  the  coronal  hole. 

^ . 4 .  The  Spectra  of  the  Second  Pe  r lod  Corona  1  Holes 

Determining  a  spectrum  for  a  coronal  hole  Is  a  rather 
difficult  task  when  there  are  many  active  regions  on  the 
solar  disk.  The  spectrum  Is  actually  the  ratio  of  the 
brightness  temperature  Inside  the  coronal  hole  at  a 
particular  wavelength  to  the  brightness  temperature  that  we 
would  have  observed  In  the  same  region  if  the  coronal  hole 
had  not  been  present.  The  best  approach  for  determining  the 
values  for  the  quiet  sun  is  to  measure  the  brightness 
temperature  In  a  quiet  area  at  the  same  distance  from  the 
center  of  the  solar  disk  and  at  the  same  latitude  as  the 
coronal  hole.  Unfortunately,  when  there  are  many  active 
regions  present,  such  quiet  areas  become  difficult  to  find, 
especially  since  the  relatively  large  beamwldths  of  the 
antennas  increase  the  area  of  the  solar  disk  affected  by  the 
active  regions. 

This  problem  was  especially  acute  during  the  second 
period,  when  two  coronal  holes,  one  In  the  northern  and  one 
in  the  southern  hemisphere,  were  situated  near  the  central 
meridian  at  approximately  the  same  distance  from  the  center 
of  the  solar  disk.  In  addition,  four  active  regions  were 
nearly  symmetrically  arranged  In  the  four  quadrants  of  the 
solar  disk.  This  made  It  unusually  difficult  to  determine 
the  appropriate  quiet  solar  brightness  temperatures,  which 
introduced  some  additional  uncertainty  In  the  determination 


(if  the  radio  spectra. 

In  L  li  i  s  section  we  will  determine  the  spectra  of  the  two 
coronal  iioles,  the  one  in  the  south  and  the  one  in  the 
north,  on  two  consecutive  days,  July  Ifi  and  17.  The 
determination  of  the  brightness  temperature  ratio  will  be 
discussed  for  eacli  map.  The  Toyokawa  naps  will  not  be  used 
because  severe  sidelobe  problems  make  them  completely 
unreliable . 

r  ii  e  Spectrum  o  i  the  Northern  Coronal  Hole  2.  Jn  1 V  1 6 »  1978 

Tlie  coronal  iiole  in  tlie  nortli  was  an  isolated  hole,  not 
connected  to  the  polar  coronal  hole.  Since  this  hole  does 
not  cover  a  wide  range  of  latitudes,  ttiere  is  no  point  in 
trying  to  get  a  spectrum  at  several  different  points  within 
tile  hole.  Indeed,  such  an  attempt  would  probably  simply 
result  in  misinformation  because  near  the  boundary  the 
antenna  beams  would  be  seeing  as  much  area  outside  the  area 
of  tlie  coronal  hole  as  inside.  We  will  therefore  look  for 
the  [loint  of  maximum  contrast. 

In  determining  the  spectrum  of  a  coronal  hole  when  we 
liad  two  maps  at  very  similar  wavelengths  we  used  the  map 
with  tlie  lilgher  resolution.  Thus,  wc  liave  chosen  the 
Haystack  map  at  1.93  cm  in  preference  to  the  La  Posta  2.0  cm 
map  and  tiie  Haystack  maps  at  1.38  cm  and  0.69  cm  in 
preference  to  tlie  Itapetlnga  maps  at  1.39  cm  and  0.70  cm. 

Table  0  lists  the  observatories;  wavelengths;  the 
quiet  sun  brightness  temperatures,  T^  in  relative  units;  the 
coronal  hole  brightness  temperatures,  Tj^  in  the  same  units; 


and  the  ratios  of  these  two  brightness  temperatures  for  the 


northern 

co-'onal  hole 

on  July  lb,  1973 

.  The 

Culgoora  map  Is 

difficult 

to  interpret 

and  the 

brightness  temperature  ratio 

at  375  cm 

must  be  treated  with 

cau  t 1  on 

.  The 

values  h  t  a  i  p  e  d 

most  likely  represent 

a  lower 

bound,  and  the 

s  pc  c  t  ra 1  value 

at  375  cm 

could  be  considerably 

closer 

to  1  . 

TABLE 

6 

Spectrum  of  the 

Northern 

Coronal 

Hole  - 

July  16 

Observatory 

A(cm) 

T 

q 

lil 

Culgoora 

375 

356 

3  06 

0.86 

Arecibo 

2  1  .4 

80 

55 

0.69 

Arecibo 

1  1  .4 

85 

7  0 

0.8  2 

Haystack 

3.8  1 

100 

92 

0.9  2 

Haystack 
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Of  particular  Interest  are  the  Haystack  observations. 
The  region  in  the  northern  hemisphere  corresponding  most 
closely  to  the  coronal  hole,  as  determined  from  ttie  I  OH  10  S 
Helium  maps  and  the  Arecibo  radio  maps,  is  a  region  of 
distinctly  lower  brightness  than  the  quiet  sun  level. 
However,  at  the  northeastern  boundary  of  the  coronal  hole 
region  the  3.81  cm  map  shows  the  suggestion  of  an 
enhancement.  Furthermore,  at  1.93  cm  and  1.33  cm  these 
regions  appear  as  very  distinct  enlianccme  nt  s  approximately 
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2-1  1)  r  i '>,li  t  f  r  tli.in  the  expected  quiet  sun  level.  This 
e  nil  .1  tice  ne  11 1  is  consistent  with  that  previously  reported  hy 
Wefei',  ct  nl.  (1976)  and  Wefer  and  Papap,  iannis  (1977)  for 
coronal  holes  around  2  cm.  Since  the  spatial  resolution  of 
the  La  Posta  telescope  used  by  Wefer  et  al.  is  much  poorer 
than  the  'lavstack  telescope,  it  is  conceivable  that  what  was 
heinp  observed  were  similar  small  enhancements  near  the 
coronal  hole  boundaries.  Although  no  plape  area  can  be  seen 
in  the  neiphhorhood  of  this  enhancement,  it  is  quite 
possible  that  near  the  boundaries  of  a  coronal  hole,  closed 
"lap  net  ic  field  lines  lead  to  higher  densities  and  hence  in 
some  cases  to  regions  of  slightly  enhanced  emission. 

r  he  S  p  e  c  t  r  u  m  o_^  t  he  Mo  r  t  he  r  n  Corona  1  Hole  2.  d  o  ^  y  1  7  ,  1  9  7  H 

'In  duly  17,  wo  have  an  additional  map  from  Culgoora, 
Australia,  hut  we  are  missing  the  0.69  cm  map  from  the 
daystack  observatory.  The  Culgoora  169  MHz  map  is 
complicated  by  sldelobe  problems,  which  make  it  impossible 
to  locate  the  coronal  holes  unambiguously  and  to  determine  a 
relative  brightness  temperature  for  them. 

Table  7  presents  the  results  for  July  17.  The  spectrum 
is  quite  consistent  with  the  results  of  July  16.  The 
contrast  appears  to  be  slightly  greater  at  11.5,  3.31,  and 
1.33  cm,  but  the  error  bars  for  these  measurements  are  about 
2-3%  so  these  differences  are  not  serious.  We  note  again 
that  the  coronal  hole  appears  as  a  depression  at  1,93  cm  and 
again  there  are  two  small  enhancements  (about  2-3%)  near  the 
northeastern  boundary  of  the  coronal  hole. 
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figure  45  shows  the  spectra  for  the  northern  coronal 
hole  for  both  days.  The  solid  line  is  the  spectrum  for  July 
16,  while  the  dashed  line  is  the  spectrum  for  July  17.  The 
error  bars  were  estimated  directly  from  the  maps. 


85 


The  Spectra  of  the  Southern  Coronal  Hole  ^  July  1 6 

The  southern  coronal  hole  consists  of  two  parts.  The 
primary  hole  is  an  extension  of  the  polar  coronal  hole, 
which  reaches  a  latitude  of  approximately  30  S.  The  other 
part  is  a  small  coronal  hole  northeast  of  the  main  extension 
at  about  10*  S  latitude. 
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This  small  region  lies  very  near  a  filament,  which  makes 
interpreting  the  maps  rather  difficult.  In  order  to  avoid 
this  difficulty,  we  will  concentrate  on  the  main  coronal 
hole.  Since  this  coronal  hole  covers  a  large  range  of 
latitudes,  we  can  gain  some  additional  information  by 
determining  the  coronal  hole  spectrum  not  only  at  the 
position  of  greatest  contrast,  but  also  at  a  larger  aiming 
parameter  near  the  limb. 

Table  8  presents  the  data  for  two  different  positions 
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in  the  coronal  hole.  The  first  is  at  the  region  of  greatest 
contrast  (about  40  S)  and  the  second  is  near  the 
photosphoric  limb.  In  general,  near  the  limb  the  brightness 
contrast  is  much  less.  However,  it  is  interesting  that  the 
Haystack  map  at  0.69  cm  and  the  Itapetinga  map  at  0.70  cm 
suggest  a  slightly  increased  contrast  between  the  coronal 
tiiile  and  the  normal  brlglitness  near  the  limb.  The  Haystack 
in.i  p  It  1.91  c  ;n  shows  very  little  contrast  near  the  limb, 
except  for  two  small  regions  near  the  coronal  hole.  The  one 
on  the  left  shows  a  depression  of  about  2%,  while  the  one  on 
the  riglit  shows  an  enhancement  of  about  1-2X.  This  slightly 
enhanced  region  (whiclt  is  almost  directly  on  the  central 
leridian)  may  not  he  within  the  area  of  the  coronal  hole  and 
nay,  therefore  be  a  region  similar  to  the  sli  ’ht 
■  •  n  1  I ..  c  e -u' II  t.  s  near  the  northern  coronal  hole.  If,  however, 
this  region  is  within  ttie  coronal  hole,  the  result  would  be 
consist  ttnt  with  the  earlier  results  discussed  in  Chapter  3, 

The  Spectra  o  f  t  he  Southern  Coronal  Hole  -  July  1 7 

i'ahlc  9  gives  tlie  coronal  hole  spectra  at  the  points  of 
highest  contrast  and  near  the  limb  for  July  17.  The  Culgoora 
160  tlilx  map  is  sufficiently  clear  in  the  southern  hemisphere 
(as  opposed  to  the  norttiern  hemisphere  -  see  above)  to  allow 
us  to  include  it  in  the  determination  of  the  spectrum.  As 
with  the  northern  coronal  hole,  the  r'sults  are  similar 
enough  to  the  results  for  July  16  to  be  considered 
consistent.  Figure  46  shows  the  spectra  for  the  southern 
coronal  hole  for  both  days  and  both  positions. 
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Chapter  5 

SPKCTRA  ANR  CORONAL  MOLE  MODELS 
5.1.  The  Spectra  o f  Corona  1  Holes 

Over  the  period  covered  by  this  contract  we  liave  been 
able  to  determine  radio  spectra  for  four  different  coronal 
lioles.  In  one  case,  the  southern  coronal  hole  in  July,  1  9  7  5  , 
we  were  even  able  to  determine  the  spectrum  at  two  different 
positions  inside  the  hole,  one  near  the  polar  limb  and  the 
other  at  the  point  of  highest  contrast  at  about  40“S 
latitude. 

Although  the  individual  spectra  are  discussed  in  the 
previous  chapters,  it  is  instructive  to  compare  all  of  them 
here  to  see  if  there  are  any  significant  differences  among 
Che  different  coronal  holes.  Overall,  the  results  are  quite 
consistent,  and  the  differences  between  coronal  holes  seem 
to  be  within  the  error  bars  of  our  observations  and 
analysis.  The  greatest  area  of  disagreement  occurs  around 
2.0  cm.  The  spectrum  for  CKl  of  197  3  and  tlie  spectrum  for 
the  coronal  hole  of  September,  1977  both  show  a  slight 
enhancement  of  approximately  2%  around  ttiis  wavelength.  On 
the  other  hand,  neither  of  the  two  coronal  holes  in  the 
July,  1973  data  show  this  enhancement,  although  in  one  case 
there  was  a  region  near  and  possibly  inside  the  coronal  hole 
that  does  show  an  enhancement.  As  was  pointed  out  in  the 
previous  chapCet,  Che  Haystack  data  for  July,  1978  do  show  a 
region  of  enhanced  emission  near  the  nortliern  coronal  hole 
but  not  inside  tl>e  hole,  while  the  southern  cornal  hole 
shows  a  slight  enhancement  near  the  limb  which  may  or  may 


not  he  inside  the  coron.sl  hole. 

It  is  i  n  t  e  r(' s  l  I  n,i>  to  note  that  tlie  data  for  Clll  indicate 
a  depression  of  about  3%  at  2.tt  cm  (Furst  and  Mirth,  1975) 
and  a  depression  of  about  2%  at  3.8  cm.  This  rather  odd  dip 
in  tlie  spectrum  is  suspicious  and  may  indicate  that  the 
errors  In  determining  tlte  relative  brightness  temperatures 
were  , greater  than  estimated.  If  so,  the  data  are  consistent 
with  no  detectable  difference  in  brightness  temperature 
around  2  cm  between  a  coronal  hole  and  the  quiet  solar 
h  .a  c  k  p,  round  . 

Finally,  it  must  be  pointed  out  that  the  spectrum  near 
the  limb  for  the  July,  1978  southern  coronal  hole  shows  much 
less  contrast  wtien  compared  with  the  spectrum  of  the  same 
coronal  iiole  on  the  disk.  This,  however,  is  to  be  expected, 
since  near  tlie  limb  we  do  not  reach  as  deep  into  the  solar 
atmosphere,  whore  the  effects  are  more  pronounced,  as  we  do 
near  the  center  of  the  solar  disk. 

5.2.  I o '! o  1 1  n g  -  Procedures 

'laving  determined  the  radio  spectrum  of  a  given  coronal 
hole,  we  now  want  to  use  it  to  obtain  the  changes  that  occur 
in  the  electron  density,  temperature,  and  neutral  particle 
density  profiles  Inside  the  coronal  hole.  The  brightness 
temperature  at  any  given  wavelength  Is  a  function  of  the 
temperatures,  electron  densities,  and  neutral  particle 
densities  along  the  ray  path  to  tnc  observer.  At  shorter 
wavelengths  the  emission  originates  deeper  inside  the  solar 
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aLmobiJiieie,  while  at  longer  wavelengths  It  originates  hlglier 
up.  With  observations  available  from  0.7  cm  to  175  cm,  i.e. 
over  more  than  9  octaves,  we  are  able  to  gain  information 
from  as  deep  as  the  lower  chromosphere,  through  the 
important  transition  zone  and  all  the  way  to  the  outer 
corona . 

The  process  of  analyzing  the  data  consists  essentially 
of  using  a  ray  tracing  procedure  to  compute  the  brightness 
temperature  of  the  sun  at  a  given  wavelength  and  radial 
distance  (aiming  parameter)  from  the  center  of  the  solar 
disk.  This  is  done  first  for  tl»e  quiet  sun  and  then  for  a 
sequence  of  modified  solar  atmospheres.  The  ratios  of  the 
computed  brightness  temperatures  for  the  modified  and  Che 
quiet  solar  atmosphere  at  a  series  of  wavelengths  are  used 
to  construct  the  radio  spectrum  for  a  particular  modified 
model  of  the  solar  atmosphere.  This  Is  then  compared  with 
the  actual  radio  spectrum  of  coronal  holes  deduced  from  our 
observations  In  an  effort  to  determine  the  modified  model 
tliat  most  closely  matches  the  observed  spectral 
characteristics  of  coronal  holes. 

This  procedure  requires  a  standard  atmosphere  for  the 
luiet  sun  which  can  then  be  modified  to  find  tlie  changes 
inside  die  coronal  holes.  The  model  we  have  used  for  the 
quiet  sun  atmosphere  is  the  Boston  University  Solar  .'lodcl 
(Papagiannis  and  Kogut;  1975,  1976)  which  follows  the 
Vernazza,  et  al.  (1973)  model  for  the  chromosphere  up  to  a 
height  of  2400  km  above  the  photosphere,  and  Is  then 
connected  using  the  method  described  by  Withbroe  (1970)  to  a 
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Dupree-Goldberg  (1967)  model  for  the  transition  zone  and  the 
corona.  The  electron  density,  ;Jg,  the  temperature,  T,  and 
the  collision  frequency,  .  of  the  model  near  tlie  critical 
transition  zone  are  shown  in  Figure  47,  while  the  complete 
listing  of  the  different  parameters  of  the  model  is  given  in 
Appendix  1. 

Tlie  brightness  temperature,  ,  is  given  by  the  integral 
along  the  ray  patii 


path 
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where  T  is  the  plasma  temperature  and  T  is  the  optical 
depth.  Since  this  is  actually  calculated  in  discrete  steps, 
the  brightness  temperature  is  approximated  by  the  expression 
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Ttie  equation  for  a  ray  trajectory  is  (Pawsey  and  Smerd; 
1  9  5  3) 
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where  r  and  0  are  sun  centered  polar  co-ordinates,  a  is  the 
.1  in  illy;  par  a  net  or  (i.e.  the  radial  distance  of  the 
otisorvation  point  from  the  center  of  the  solar  disk),  and  p 
is  the  real  part  of  the  complex  index  of  refraction.  Both  r 
and  a  arc  measured  in  terms  of  the  photosphoric  radius,  11^. 
The  optical  depth  is  related  to  the  absorption  coefficient, 
K,  by  the  expression 


At  =  (c  as 

m  mm 


m  m 
cos  e 


wiiere  As  is  the  path  length  through  layer  m. 

m 

The  absorption  coefficient, k  ,  in  layer  ra,  is  a  function 

m 

of  tlie  plasma  frequency  ,  ^,m  collision  frequency, 

y  ,  in  this  layer,  and  the  radio  frequency,  f,  of  the 

c,m 

observations.  IVtien  the  imaginary  component  of  tlie  refractive 
index  is  mucii  less  than  1  and  the  real  component,  p,  is 
close  to  1  ,  tlie  absorption  coefficient  can  be  approximated 
by  the  expression 


The  coliision  frequency  is  the  sum  of  the  collision 

frequency  of  electrons  with  neutral  particles, v  and  the 

n  ,m 

collision  frequency  of  electrons  with  ions,v.  : 
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where  density  of  neutral  particles,  ‘^6,111^'’® 
electron  density,  which  is  assumnied  to  be  essentially  equal 
to  the  Ion  density,  and  T  Is  the  pla8:na  temperature.  Mtcolet 
(  1  97  3)  gives  for  the  value  »  5.4  x  10“'®  ,  wlille 
Ginzburg  (1961)  gives  an  expression  for  in  the  range  of 
3  0-50  with  only  a  logarithmic  dependence  on  N'g  and  T.  In 
general  the  absorption  coefficient  also  depends  In  a  tilghly 
complex  way  on  the  magnetic  field,  but  for  the  frequencies 
of  interest  here  and  the  magnetic  field  strengths 
characteristic  of  coronal  holes,  the  magnetic  field  effects 
are  essentially  negligible,  and  can  therefore  be  neglected. 

Once  the  temperature  and  particle  density  profiles  are 
known,  the  ray  tracing  procedure  (equation  4)  can  be  used  to 
follow  the  path  of  the  ray  and  the  brightness  temperature  of 
the  emerging  ray  will  be  given  by  equation  2.  In  order  to 
account  for  different  changes  of  the  quiet  atmospliere  model 
In  different  height  ranges  inside  the  cornal  hole,  we  have 
divided  the  standard  atmospheric  model  Into  five  regions: 
the  photosphere  (0  -  600  km),  the  Inner  chromosphere 
(600-2300  km),  the  outer  chromosphere  (2300  -  2465  km),  the 
transition  zone  (2465  -  2520  km),  and  the  corona  (2520  - 
700,000  km).  Since  our  radio  data  provide  no  information 
about  the  photosphere,  this  region  played  no  role  in  our 
analysis  and  therefore  was  left  unmodified.  The  upper  four 
regions  were  given  the  numerical  designations:  lower 
chromosphere  (0),  upper  chromosphere  (1),  transition  zone 
(2),  and  corona  (3). 

In  each  of  these  upper  four  regions,  the  electron 


density  of  each 
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density  can  be  i;iodified  by  multiplying  the 
layer  in  a  region  by  a  numerical  factor  (l»0,l,2,3). 

Similarly,  the  temperature  of  each  region  can  be  modified  by 
multiplying  the  temperatures  by  a  factor  6^  .  For  the 
modification  of  the  neutral  particle  density  we  use  a  single 
factor,  Y«  for  regions  because  their  effects  are 

important  at  only  the  lower  regions  which  are  reached  only 
by  the  shorter  wavelengths.  A  different  form  of  modification 
is  achieved  by  ctianging  the  thickness  of  each  of  the  upper 
four  regions.  The  thicknesses  are  clianged  by  multiplying 
tliem  with  the  factors  (1«  0  -  3),  which  in  effect  change 

the  gradients  of  all  quantities  in  each  region.  If  tiie 
modifications  in  two  adjacent  regions  are  different,  the 
corresponding  physical  quantities  are  interpolated  near  the 
boundary  so  as  to  assure  a  smooth  transition  from  one  region 
to  the  next.  When  all  the  numerical  factors  mentioned  above 
arc  equal  to  unity  we  simply  have  the  standard  solar  model. 

The  entire  procedure  of  modifying  the  standard 
atmospheric  mo<iel,  tracing  ray  paths  and  obtaining 
bri'jhtness  temperatures  at  several  wavelengths,  normalizing 
the  brightness  CeoperaCures  to  the  corresponding  values  for 
tlip  unmodified  atmosphere,  and  finally  plotting  the  entire 
radio  spectrum  for  a  particular  modification  of  the  standard 
model  is  accomplished  by  a  single  computer  program.  This 
program  is  listed  in  Appendix  II  along  with  a  sample  Input 
and  output. 

Although  It  Is  possible  to  determine  the  effect  of 
changing  a  single  parameter,  in  general  these  parameters 


Interact  wltU  each  otfier  in  non-lIncar  ways,  .is  is  ilie  c.ise 
between  electron  density  and  temperature.  \  particnl.ar 
example  is  tlio  relation  iictween  tlie  electron  density  and  the 
depth  of  ttie  transition  zone.  The  total  nii'iber  of  elect  ronr; 
In  the  transition  zone  does  not  ,i  .ape.ar  to  e.i.in  .•  v.- r  . 
s  i  n  I  f  i  c  1  nt  I  V  as  it  expands  in  .1  coron.il  'i>le.  I'luis  wlien  '  In- 
transition  zone  oxp.inds  l>y  a  factor  of  about  ')  or  ii  ,  as  w.is 
sue  jested  by  XlIV  data  from  Sky  lab  (Roblin,  et  al.,  1 7 ; 
duber,  et  al.,  1974;  Wlthb  roe,  1  977  ),  the  electron  nuinber 
density  decreases  to  about  20%  of  the  norm.il  value.  Tlie 
effect  of  the  different  modi  f  Icat  1  ons  on  the  ob.se  rvcd 
briahtnoss  depends  also  on  the  ainlny,  parameter,  a,  since  .it 
different  values  of  the  aiming  parameter  the  rays  follow 
different  paths  through  the  solar  atmosphere,  dear  the  lii.ib, 
for  example,  at  longer  wavelengths  we  have  the  phc n ome no n  of 
limb  brightening,  which  tends  to  obscure  the  effects  of 
coronal  holes. 

It  must  be  mentioned  that  the  entire  modeling  procedure 
..e  have  followed  assumes  a  spherically  symmetric  solar 
atmosphere  and  neglects  changes  toward  the  polar  regions, 
wiilcli  are  actually  nresent.  We  are  mostly  interested, 
however,  in  the  modifications  of  a  given  profile  of  the 
solar  atmosphere  inside  a  coronal  hole.  The  use  of  a 
somewhat  different  profile  In  the  polar  regions  would  simply 
produce  so.newhat  different  .mod  i  f  1  ca  1 1  ons  more  appropriate 
for  this  model,  but  overall  the  changes  are  likely  to  be  in 


the 


s  a  me 


direction  as  in  the  case  of  the  st.inilard  lode  1 


9H 


)  .  J  .  Mod  e  1  i  lit;  -  1< c  s  II  L  t  s 

i’o  obtain  a  model  for  the  solar  atmosphere  inside  a 
coronal  hole,  we  begin  by  determining  the  effects  on  the 
spectrum  ot  changing  each  parameter  of  the  standard  model 
separately.  These  effects,  of  course,  will  be  somewhat 
different  for  different  values  of  tiie  tlie  aiming  parameter, 
a.  If  chan,',ing  a  single  parameter,  such  as  or  T,  doesn't 
produce  l  lie  desired  results  we  try  cominations  of  the  two 
parameters  forming  .i  p.  r  i  d  of  different  numerical  values  for 
a  anti  8.  The  fitting  of  the  coronal  hole  spectrum  proceeds 
from  the  longer  wavelength  end  because  this  corresponds  to 
the  outer  .Host  layers  and  is  not  affected  by  clianges  in  tlie 
lower  re'.ions.  We  then  proceed  to  the  shorter  wavelengths 
anti  the  deeper  layers  of  tlie  atmosphere.  Once  a  reasonably 
:’,ood  fit  is  found  for  tlie  entire  spectrum,  tlie  parameters 
'•111  he  further  miidified  by  small  amounts  to  improve  tlie  fit. 

fo  demonstrate  how  this  procedure  works  we  will  consider 
i a  detail  the  process  as  it  applied  to  one  particular 
coronal  hole,  the  southern  coronal  hole  of  July  1978.  Figure 
*  1  ■>  i  o  V  s  the  computed  brightness  temperature  ratios  and  the 
observe  1 1  ratios  for  the  southern  coronal  hole  of  duly  17, 
19/.-).  fhe  pro)  i  llinary  estimates  for  the  parameters  are  given 
in  T  a  1 1 1 1’  1  . )  . 
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Preliminary  P.nramcter  Estimates 
Southern  Coronal  Hole  Disk  Spectrum  -  July  17,  1978 


a 

6 

Y 

0 

0.8 

0.9 

0.8 

2.5 

1 

0.8 

0.3 

0.8 

2.0 

2 

0.15 

0.45 

0.8 

6  .0 

3 

0.3 

0.3 

0.8 

1  .0 

In  Figures  49  -  56  we  show  the  procedure  used  to  refine 
the  original  estimates.  Each  graph  contains  4  curves  plus 
the  measured  spectrum.  In  each  case,  two  parameters  in  one 
region  (e.g.  and  )  have  been  modified  around  the 

Initial  estimates.  From  this  set  of  graphs  the  effects  of 
varying  the  different  parameters  around  the  initial 
estimutes  can  be  iletermined.  The  dominant  effect  of  changing 
parameters  In  the  chromosphere  (regions  0  and  1)  is  observed 
at  wavelengths  sliorter  than  about  5  cm.  Varying  the 
parameters  in  the  transition  zone  (region  2)  primarily 
affects  the  wavclengtti  range  from  about  5  cm  to  50  cm,  while 
chan'gin,;  the  corona  (region  3)  affects  primarily  the 
wavelengths  lonj'er  than  50  cm  but  of  course  also  has  some 
affect  on  wavelengths  in  the  5  cm  to  50  cm  range. 

One  significant  effect  is  that  of  changing  the  thickness 
of  a  region  (parameter  r,  )•  Increasing  the  thickness  of  a 
rev,  ion  wt\ilo  keeping  tl»e  density  and  temperature  unchanged 
results  in  an  enhanced  brightness  temperatui'e  at  wavelengths 
wlilch  emanate  primarily  from  the  modified  region.  However, 
if  at  the  same  time  the  density  is  decreased  so  as  to 
maintain  a  constant  number  of  particles  in  the  region,  the 
two  effects  tend  to  cancel  out,  and  only  a  s ma 11  change  in 
brlgiitness  results. 
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observed  spectrum 


After  observing  the*  effects  produced  by  cban;»,  ing 
different  pairs  of  parameters,  improved  fits  to  tlie  observed 
spectrum  can  be  made  by  varying  all  tbe  parai:ieters  over  a 
limited  range.  The  result  of  this  final  refinement  of  the 
liarameters  is  shown  in  Figure  S7  and  the  values  of  the 
parameters  are  given  in  Table  11. 


TABLF.  11 

Final  Parameter  Estimates 

Southern  Coronal  Hole  Disk  Spectrum  -  July  17,  1978 
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Our  best  results  for  the  coronal  hole  modeling  come  from 
the  July,  1973  observations.  We  have  used  these  data  to 
produce  models  of  both  the  southern  coronal  hole  and  the 
nortiiern  coronal  hole  for  both  July  16  and  July  17.  In 
addition,  we  have  produced  separate  models  for  the  two 
different  aiming  parameters  for  the  southern  coronal  hole. 
These  results  are  summarir.ed  in  Figures  57  -  62.  In  eacii 
figure  the  observed  spectrum  Is  shown  as  open  boxes,  while 
the  theoretical  spectrum  is  shown  as  open  circles.  As  can  be 
seen,  the  models  fit  the  observations  quite  well.  Tlie  values 
of  the  parameters  for  eacli  model  are  listed  in  Table  12.  We 
note  that  the  parameters  are  quite  consistent  both  for  the 
two  different  days  and  for  the  two  coronal  lioles. 
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The  layers 

of  the  solar  atmosphere 

which 

c  on  t 

r  1  b  u  t  e 

the 

l.iri'est  fr.iction  of  the  radio 

flux  observed, 

vary 

f  ro'.n 

w  a  v  e  1  c  n  t  h  to 

wavelength.  Figures  63 

a  nd 

64 

show 

the 

c  on  t  r  1  bvi  L  t  on 

to  the  flux  as  a 

funct Ion 

of  helglit 

for 

different  wavelengths.  These  fisurcs  are  derived  from  the 
•nodel  of  t!ie  stuithern  coronal  hole,  but  usinv,  a  slightly 


different  modification  for  a  different  coronal  hole  does  not 
si", nlflcantly  alter  the  contribution  profiles. 


SOUTHERN  CORONRL  HOLE  DISK  SPECTRUM  -  JULY  17 


FIGURE  57 


CORONRL  HOLE  SPECTRUM  -  JULY  17 
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SOUTHERN  CORONAL  HOLE  LIMB  SPECTRUM  -  JULY  17 


SOUTHERN  CORONAL  HOLE  DISK  SPECTRUM  -  JULY  16 


NORTHERN  CQR0NRL  HOLE  SPECTRUM  -  JULY  16 


FIGURE  62 


CONTRIBUTION 
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AIMING  PARAMETER -0.70 


FIGURE  63 


CONTRIBUTION 
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AIMING  PARAMETER  »0.90 


FIGURE  f>k 


5.4 


Conclusions 


The  models  we  have  computed  for  coronal  holes  are  cjiiite 
consistent  from  one  hole  to  another  and  from  one  viewing 
ani^le  (aiming  parameter)  to  another.  The  general  features 
are: 

1)  The  chromosphere  expands  in  thickness  hy  about  a 

factor  of  2,  but  apparently  fills  In  with  material 
from  the  photosphere,  since  the  values  for  and 
Ng  both  decrease  by  only  about  '2.0Z  .  The 
controversial  enhancement  near  2  cm  can  be 
obtained  through  a  significant  increase  of  the 
neutral  particle  density  in  this  region. 

2)  The  transition  zone  expands  In  thickness  by  about  a 

factor  of  5  to  6  and  the  electron  density  drops  by 
a  similar  factor. 

3)  The  corona  Itself  does  not  appear  to  expand,  but  is 

depleted  in  particle  density. 

4)  The  temperatures  are  typically  reduced  by  about 
20%,  except  In  ttie  transition  zone,  where  the 
reduction  is  around  50%.  This  indicates  continued 
heating  of  these  regions  despite  their  expansion. 

The  coronal  hole  models  obtained,  together  with  other 
relevent  information,  such  as  solar  wind  speeds,  etc.,  can 
now  be  used  to  develop  a  quantitative  understanding  of  the 
physical  processes  operating  inside  coronal  holes  and  to 
bring  theory  and  observations  into  full  agreement.  This  Is  a 
task  that  was  beyond  the  goals  of  this  project,  but  we  hope 
to  be  able  to  undertake  It  In  the  near  future. 
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APPENDIX  1 

STANDARD  SOLAR  ATMOSPHERE  MODEL 


HT(km) 

TCK) 

N(el/cm’) 

N« (par t/cm’) 

487720.0 

1718400.0 

0.76684e 

07 

O.OOOOOe 

00 

457720.0 

1759800.0 

0.95640e 

07 

O.OOOOOe 

00 

587720.0 

1762500,0 

0.11620c 

08 

O.OOOOOe 

00 

557720.0 

1776500.0 

0.15341e 

08 

O.OOOOOe 

00 

527720.0 

1791200,0 

0.15426e 

08 

O.OOOOOe 

00 

297720.0 

1806500.0 

0,17975e 

08 

O.OOOOOe 

00 

267720.0 

1822400.0 

0.21118e 

08 

O.OOOOOe 

00 

257720.0 

1859000  .0 

0,25037e 

08 

O.OOOOOe 

00 

207720.0 

1856200.0 

0.29977e 

08 

O.OOOOOe 

00 

177720.0 

1874200.0 

0,56282e 

08 

O.OOOOOe 

00 

147720.0 

1895000,0 

0, 44434e 

08 

O.OOOOOe 

00 

117720.0 

1912700.0 

0. 551306 

08 

O.OOOOOe 

00 

87720.0 

1955500.0 

0.69391e 

08 

O.OOOOOe 

00 

77720.0 

1940400.0 

0, 751816 

08 

O.OOOOOe 

00 

67720.0 

1947600.0 

0.816056 

08 

O.OOOOOe 

00 

57720.0 

1955000.0 

0,887436 

08 

O.OOOOOe 

00 

47720.0 

1962500,0 

0.96696e 

08 

O.OOOOOe 

00 

57720.0 

1970100.0 

0. 150006 

09 

O.OOOOOe 

00 

27720.0 

1959700.0 

0. 153096 

09 

O.OOOOOe 

00 

17720.0 

1692900.0 

0. 177216 

09 

O.OOOOOe 

00 

16720.0 

1659900.0 

0. 180746 

09 

O.OOOOOe 

00 

15720.0 

1625200,0 

0,184606 

09 

O.OOOOOe 

00 

14720.0 

1588600,0 

0. 188856 

09 

O.OOOOOe 

00 

15720,0 

1549700.0 

0.195586 

09 

O.OOOOOe 

00 

12720.0 

1508400.0 

0.198896 

09 

O.OOOOOe 

00 

11720,0 

1464100.0 

0, 204916 

09 

O.OOOOOe 

00 

10720.0 

1416200,0 

0.211836 

09 

O.OOOOOe 

00 

9720.0 

1564000.0 

0,219946 

09 

O.OOOOOe 

00 

8720,0 

1506400,0 

0. 229646 

09 

O.OOOOOe 

00 

7720,0 

1241800,0 

0.241586 

09 

O.OOOOOe 

00 

6720.0 

1167700,0 

0.256916 

09 

O.OOOOOe 

00 

5720,0 

1080000,0 

0,277796 

09 

O.OOOOOe 

00 

4720.0 

970070,0 

0.309266 

09 

O.OOOOOe 

00 

5720,0 

816610,0 

0,367376 

09 

O.OOOOOe 

00 

2720.0 

505800,0 

0.595486 

09 

O.OOOOOe 

00 

2620.0 

428270.0 

0,700496 

09 

O.OOOOOe 

00 

2520.0 

286250.0 

0, 104806 

10 

O.OOOOOe 

00 

2510.0 

257990,0 

0.116286 

10 

O.OOOOOe 

00 

2500.0 

219070,0 

0.156946 

10 

O.OOOOOe 

00 

2495.0 

190820,0 

0,157226 

10 

O.OOOOOe 

00 

2490.0 

145250,0 

0. 206546 

10 

O.OOOOOe 

00 

2489.0 

150290,0 

0.250256 

10 

O.OOOOOe 

00 

2488.0 

109290.0 

0,274506 

10 

O.OOOOOe 

00 

2487.0 

65000.0 

0,694206 

10 

0,999426 

06 

2485.0 

45000.0 

0.970006 

10 

0. 599656 

07 

2484.0 

54000,0 

0. 125906 

11 

0. 399976 

08 

2485.0 

27000,0 

0,149206 

11 

0, 210006 

09 

2482.0 

24000.0 

0, 157706 

11 

0.520006 

09 

2481.0 

22500,0 

0. 157906 

11 

0,830006 

09 
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HT(km) 

ICK) 

N{e1 /cm’ ) 

Nr,(part/cm’  ) 

2480.0 

22000.0 

0.15820e 

11 

0. 960006 

09 

2465.0 

21000.0 

0.16420e 

11 

0,102006 

10 

2445.0 

20600,0 

0.17090e 

11 

0, 970006 

09 

2415.0 

20400,0 

0. 178306 

11 

0,940006 

09 

2375.0 

20300.0 

0.18660e 

11 

0. 970006 

09 

2315.0 

20200.0 

0. 196706 

11 

0,116006 

10 

2305.0 

19800.0 

0,201106 

11 

0, 139006 

10 

2302.0 

16000,0 

0, 226706 

11 

0,337006 

10 

2298.0 

13000.0 

0. 250206 

11 

0, 636006 

10 

2295.0 

8800.0 

0. 289906 

11 

0.139306 

11 

2282.0 

8550.0 

0.28360e 

11 

0. 177606 

11 

2268.0 

8300.0 

0. 276006 

11 

0.218806 

11 

2235.0 

7800,0 

0.26260e 

11 

0.322306 

11 

2170.0 

7300,0 

0. 245906 

11 

0,528806 

11 

2095.0 

6940,0 

0,229706 

11 

0.841606 

11 

2030.0 

6720.0 

0. 218406 

11 

0. 122076 

12 

1945.0 

6620.0 

0.222106 

11 

0.185546 

12 

1810.0 

6560.0 

0. 252806 

11 

0.358906 

12 

1690.0 

6500.0 

0,293206 

11 

0.665086 

12 

1560.0 

6440,0 

0, 374806 

11 

0.135266 

13 

1430.0 

6390.0 

0,515006 

11 

0.287746 

13 

1280,0 

6330.0 

0.793806 

11 

0,726436 

13 

1170.0 

6230.0 

0,105506 

12 

0,148506 

14 

1030,0 

5930.0 

0.12690e 

12 

0.398856 

14 

980.0 

5760.0 

0,121806 

12 

0. 586946 

14 

930.0 

5560,0 

0.11060c 

12 

0,881706 

14 

880.0 

5360.0 

0. 972806 

11 

0.134136 

15 

830.0 

5150.0 

0, 855606 

11 

0.207156 

15 

780.0 

4890.0 

0, 779106 

11 

0.329976 

15 

730.0 

4600.0 

0, 842406 

11 

0. 542696 

15 

680.0 

4350.0 

0.112206 

12 

0.913506 

15 

630.0 

4170.0 

0,161506 

12 

0.155606 

16 

605.0 

4150.0 

0,196106 

12 

0. 200606 

16 

580. C 

4150.0 

0,241006 

12 

0. 257606 

16 

530.0 

4200.0 

0. 375406 

12 

0.418706 

16 

480.0 

4350.0 

0. 615106 

12 

0. 657306 

16 

430,0 

4500.0 

0, 990406 

12 

0. 101606 

17 

380,0 

4640.0 

0. 156206 

13 

0. 155406 

17 

330.0 

4780.0 

0.243106 

13 

0. 234606 

17 

280.0 

4930.0 

0, 376106 

13 

0. 349406 

17 

230.0 

5100.0 

0. 584606 

13 

0.511706 

17 

180.0 

5370,0 

0,964506 

13 

0. 724396 

17 

130.0 

5750.0 

0,189806 

14 

0,985746 

17 

100.0 

6070.0 

0.348506 

14 

0.115486 

18 

80.0 

6300,0 

0.551106 

14 

0,127366 

18 

65.0 

6600,0 

0. 954906 

14 

0.134036 

18 

50.0 

7100,0 

0, 220906 

15 

0.136416 

18 

30.0 

8000,0 

0, 804006 

15 

0.134336 

18 

0,0 

9000,0 

0,260106 

16 

0.134956 

18 
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Al’PKNDlX  II 
RAY  TRACING  PROGRAM 


/ 1  oad 
/  pa  rm 
/job 
c 

cccc 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


1  4 
7  1 


c 

c 

cccc 

c 

c 

c 

c 


fortgl ,nsegs=6 
deck 
nogo 


program 
dime  ns  Ion 


" s  o 1 mf  w" 


dh(200)  ,h(2()0)  ,  t(200)  ,cr(2o0)  ,dtau(200)  ,crm(2l)  , 
ar(20n),br(200),dr(200),bt(2l),btq(21),spec(21) 
r,rmral,s,com(40),rcrin(21),ri!i,ocr(51),rl,r2,r3,rlpl, 
r2pl  ,r3pl,p,q,rmpl  ,ppl  ,pp2,rnl  ,rt,r3nl  ,r4,r4nl ,r0, 
rOpl ,rlm3 

1(  2  I  )  ,  nu(2 1 ) , nn(200) ,mn2(  2  00) , Inu, ne(200) , htel (2 1 ) . 
hcrl(21),hcr2(21),nue(200),nuep(20()),nuen(200), 
tq(  200) ,ncq( 2  00) , ro2(2  00),hr(21),tauref(2l),the(200), 
theinf(2l),er(2  00),xr(20  0),nnq(2  00),liq(20  0) 
logical  qflag/. false. /,aflg/. false./ 
lnteger*A  gtltl(15) 


Intege  r 


real 


r m» 1 8  1 

rm  «  no.  of  levels  In  tbe  solar  atmosphere, 
r  mp  1  =•  rm+  1 
rmml -rm- 1 
xc-2 . 9  9  79e  1 0 

xc-  speed  of  light  (cm/s), 
r sun-6 . 9  59  9e  5 

rsun  »  radius  of  the  sun  (km). 
pl-3. 141 5927 
pi  -  that's  right,  pi. 
crd-180.0/pl 

crd  -  conversion  of  radians  to  degrees. 
x5- 1 .0e5 

x5  -  number  of  cm/km. 
xk  f c- 1 . 1 e-6 

xkfc  -  constant  for  computing  conductive  flux. 
xh3-3 . ia260e9 

xb3  -  constant  for  computing  the  plasma  frequency. 

do  14  1-1,51 
oc  r ( 1 ) -0 
do  71  1-4,36,4 
ocr(l+l )-i/4 
ocr(41)-l 
x9-l .Oe-9 


input  wavelengths,  almpar,  spectrum,  and  region  limits  cccccccccc 
read  ( 5 , * )  Im 


Im  is  the  number  of  wavelengths  that  are 


to  be  processed 


\ 


4 


do  16  1-1,  im 
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16  read  (5,617)  l(i) 

617  format  (3fl0.0) 

c  l(i)  “  wavelength  (cm). 

c  Im  °  number  of  w.ivelengths  to  process. 

HOI  read  (5,602,end“998)  almpar,spec 

602  format  ( f 7 . 0 , 2  I f 3 . 2 ) 

read  (5,851)gtltl 
851  format(15a4) 

aflg~. false. 

c  spec(i)  “  coronal  hole  spectrum. 

d2«ain\par*airapar 

read  (5,611)  rO,rl,r2,r3,lp,alfaO,alfal,alfa2,alfa3,betaO, 

1  betal , beta2 , beta! , gamma 

611  format(5i3,9f6.0) 

c 

c 

if(qflag)go  to  802 

read  (5,601)  (coir(l),l«l,40) 

qfla>»=.true. 

6Ul  format  (40a2) 

c 

do  1  r “ 1  , r  m 

read  (  5  ,  6 00  , e nd  =  993  )  hq ( r ) , t q ( r ) , ncq ( r ) , n nq ( r ) 

600  format  ( f 8 . 0 , f 1 0 . 0 , e 1 2 .0 , e 1  1  . 0 ) 

c  l)q(r)  »  altitude  of  layer  r  (km), 

c  tq(r)  =*  temperature  in  layer  r  (k). 

c  ncq(r)»  it  electrons  per  cm**3 

c  nnq(r)“  9  neutral  atoms  per  cm**3. 

1  continue 

c 

802  alfa0»1.0 

alfal»l  .0 
a  I  f  a2 “ 1 .0 
alfa3»l  .0 
b  e  t  a  0  » I  .0 
betal » I  .0 
b  e  t  a  2  «  1 .0 
b  e  t  a  3  =  1  .0 
ze  ta0»  I  .0 
ze  ta  1  =■  1  .0 
zeta2»‘l  .0 
ze  ta3» 1  .0 
gamma>l  .0 
lcount"0 

the(rmpl)»crd*arsln(alrapar*rsun/(h(rmpl)+rsun)) 

c 

cccc  modify  the  quiet  sun  model  c c cccccccccccc cccc cccccc cc cccc cccccc c c 
c 

7  rOpl-rO+1 
r  1  pl«rl  1-1 
r2pl-r2+l 
r  3  pi  •  r 3+ 1 
r3ml*r3-  1 
r4“ r 3- 2 


-  129- 


r Am  1 »  r 4- 1 
r  1  m  3*  r  1  -  3 
dr>  9  r»l,rO 
t ( r )»tq ( r  ) 

9  ne(r)“neq(r) 
c 

t  p»  t  q  (  rO  ) 

do  6  r=rOpl , r 1 

t(r)»tp+betat)*(tq(r)-tp) 

6  ne(r)=neq(r)*alfaO 
c 

tp*tq( rO) 

do  8  r=rlpl,r2 

t(r)=tp+betal*(tq(r)-Cp) 

if  ( t ( r ) . Rt . t q ( r 2 ) )  t(r)»tq(r2) 

3  ne(r)“neq(r)*alfal 
c 

tp-t (r2 ) 

do  19  r*r3,r3i 

t(r)=tp+heta3*(tq(r)-tp) 

19  ne( r ) = neq ( r ) *a I f a3 
t(rmpl  )=«t(rm) 
c 

tp=t(r2) 

do  17  r*r2pl,r3inl 
t(r)=tp+beta2*(tq(r)-tp) 

17  ne ( r )= neq ( r ) *a If a2 

c 
c 

c  recompute  the  altitude  of  each  level 

c  this  also  changes  the  gradients  in  each  region 

c 

do  900  r=l , rO 

900  ti(r)*hq(r) 
hl“hq ( rO) 

do  901  r=r0pl,rl 

h(r)”hl  +  zeta0*(hq(r)-hq(r-l)  ) 

901  hl-h(r) 

do  902  r-rl pi , r2 

h(r)*hl  +  zetal*(hq(r )-hq(r-l )  ) 

902  hl-h(r) 

do  903  r»r2pl,r3ml 

h(r)»lil  +  zeta2*(hq(r)-hq(r-l)) 

903  hl-h(r) 

do  904  r*'r3,rra 

h(r)"hl+zeta3*(hq(r)-hq(r-l)) 

904  hl»h^r) 

do  905  r^l.cm 
x"(rsun+h(r))/r8un 
ro2 ( r )“x*x 

905  cont  Inue 
c 

do  2  r “2 , rmm 1 
dh(r)-(h(r+l)-h(r-l))*0.5 
dh( l)-h(2)-h( 1 ) 


2 


I  30 


c 


2'i 

c 

c 


2  9 
c 
c 


2f) 

c 


c 


c 

C  C  C  C' 

c 

c 

c 

c 


dll  (  rm  )  »li  (  r;n  )-ti  (  r  Tini  1  ) 

Jh(r)>'  shell  thickness  (km), 
li  (  r  111  p  1  )  =  h  (  r  m  )  +0 . 5  *  d  h  (  r  m  ) 

•natch  the  electron  density  at  r4»*>r3. 

s3=(ne(r3pl)-ne(r3))/(h(r3pl)-h(r3)) 

do  23  r“r2|>l,r3 

nep  =  ne(r3)  +  s3*(h(r)-h(r3)  ) 

If  ( ne( r ) .ge . nep)  go  to  23 
neC  r)  =  nep 
cont I nue 

match  the  temperature  at  r4==>r3. 
s4  =  (  t(r3)-t(rA) )/(h(r3)-h(rA)  ) 
do  29  r=r4,r3 

t( r)  =  t(r4)+s4*(h(r)-h(r4)  ) 

natch  Che  electron  density  at  r=  r  1  ni  3  =  =  >  r  1  +  1  . 

sl=(ne(ripl)-ne(rlm3))/(h(rlpl)-h(rlni3)) 

do  2fi  r=rlm3,rlpl 

iie(r)  =  ne(rlni3)+sl*(h(r)-h(rlni3)) 


do  11  r  =  1  , r  m 
g  a  m  =  V,  a  ii  m  a 

if  (gam.eq.l.OO)  go  to  11 
if  ( n nq ( r  )  .  le . 0 . 0 )  go  to  11 
a  If  ao  =  ne(r)/neq(r) 
betae»tq(r)/t(r) 
betae  Is  inverted! 

ga:nmax“hetae+(  2 .0*neq(r)/nnq(  r  )  )*(betae-alfae) 
If  ( gamma .  't . gammax)  gam  =  gammax 
nn(r)=nnq(r)*gam 


electron  collision  frequencies 


;  c  c  c  c  c  c  c  c 


nicolet,  in.  (1  9  53),  "the  collision  frequency  of  electrons  in  the 
Ionosphere",  journal  of  atmospheric  and  terrestrial  physics, 
vol.  3,  pp.  200-211. 


c  const.  II  ts  for  collisions  with  neutrals. 

21  xkn  =  2 .6022ef) 

s  1  >'  =  2 . 0  e  -  8 

c  slg  »  collision  distance  for  electrons  with  neutrals  (cm). 

-s  i  2  *  s  1  g  *  s  i  g 

c  constants  for  collisions  with  positive  ions. 

x32=(-3 .0/2 .0) 
x23-(-2 .0/3 .0) 
x  b 1 « I  .81668 
xk4-5 .7  2965e6 
c 

do  15  r ”  I  ,  r m 

nuen(r)“xkn*nn(r)*8lg2*sqrt(t(r)) 

c.  nuen(r)  -  collision  frequency  between  electrons  and  neutrals. 


al2-alQg( 1 .0  +  (xk4*t(r)*t(r)*(2.0*ne(r))**x23)  ) 


131 


c 

c 

1  5 
c 
c 

c  c  c  c 
c 


603 


c 

3  1  1 


6  04 

615 

22 

c 

c  c  c  c 
c 

c  c  c  c 


c 


5 

c 

cc 

cccc 
c  c 
c 


c 


nuep(r)»xbl*al 2*ne( r)*( t(r)**x32) 

nuep(r)  -  collision  frequency  between  electrons  and  protrons. 
nue( r)“nuep(r)+nuen(r) 

nue(r)  •  electron  collision  frequency. 

1 f ( 1 p . eq . -  1 ) go  to  311 

wrlte(ll,616)(com(n.j“lt^O),rO,alfaO,betaO,zeta(),rl,alfal, 

1  betal ,zctal ,r2,alfa2,bcta2,zeta2,r3,alfa3,beta3,zeta3. 


2  gamma, aimpar 
write  (11,603) 
f  o  rma  t  ( 

I0  5h  r 

h 

dh 

t 

ne 

1  n  n  n  11  e 

niiep 

nuen 

ratio 

f 

p*p/f /3 

2,102h  (-)  (km)  (km)  (k)  (cm**-3)  (cm**-3)  (hz) 

3(h7. )  (bz)  c/s*cm*cm  /) 

do  22  r»2,rin 
ratio=nuen( r)/nuep( r) 

fcond“xkfc*t(r)*t(r)*sqrt(t(r))*(t(r  +  l)-t(r-l))/(li(r+l)-h(r-l))‘ 
f condd«f  cond 

if  (fcondd.eq.0.00)  fcondd“1.00 

p2of”ne(r)*t(r)*ne(r)*t(r)/fcondd 

alfac=ne(r)/neq(r) 

betae*t(r)/tq(r) 

i f ( 1 p . eq . -  1 ) go  CO  22 

write  (11,604)  r , h{ r ) , dh ( r ) , t ( r ) , ne ( r ) , nn( r ) , nue( r ) , nuep( r ) , 

1  nuen(r),ratio,fcond,p2of,alfae,betae 

format  ( i 6 , f 3 . 0 , f 8 . 1 , f 9 .0 , 1 p5e 1 0 . 3  ,  5e 9 . 2  ) 

if  ( r . eq . r I . or . r . eq . r 2 . or . r . eq . r 3 . o r . r . eq . rO)  write  (11,615) 

format  (lx,13(inh . )) 

cont  inue 

optical  depth  and  contribution  function  routines  c ccc c cc ccc c c cc c c 

initialize  variables, 
do  5  i»l  ,  Im 
nu(i)*‘xc/l(l) 

nu(l)  frequency  corresponding  to  l(i)  (hz). 

bt(i)«f).0 

htel ( i  )-0 .0 

hcrl  (  1  )  =  0.() 

hcr2( 1 )»0.0 

rcrm( i  )■ I 

c  rm ( i ) “O  .  0 


wavelength  loop. 

do  999  i*l  ,  ini 
lnu“a  log(  nil  (  1  )  ) 

if  (l(i).ge.0.050)  go  to  13 


6  13 


write  (6,607)  (coni(j),j«l,A0),l(l) 

write  (6,613) 

format  (5x,  57hthe  r  ey  1  e  1  j’.h- j  e  a  ns  approximation  is  used  in  this 
lrogram./5x,  72hthis  makes  computations  at  wavelengths  less  tiian  J 
2050  cm  inappropriate.) 
go  to  999 
c 

cccc  compute  dtau  for  each  layer  and  locate  the  turning-point  level. 

13  ltl=l(l)*l(i) 

om=‘2.0*pi*nu(i) 
om2*  om*om 
mu2(rnpl)»1.00 
p  =  0 

c 

Jo  10  j  =  1  , r  m 
r  =  rmp 1  - j 
omo2=xb3*nc( r) 
yp“om2+nue(r)*nue(r) 
xmf=l .00-(omo2/yp) 
xnf=nue(  r)*om<)2  /  (om*yp) 
ypp=(xfi£*xnf  )/(xmf*xmf) 
chi2=0.25*xnf*ypp*(1.0-0.25*ypp) 
inu2(r)  =  xraf  +  chl2 
if  ( mu2 ( r  ) , 1 e . 0 . 00 )  go  to  60 
c  if  ( mu ** 2 . 1 c . 0 . 00 )  the  wave  is  reflected. 

zi=*l  .00-(d2/  (mu2(r)*ro2(r))) 
if  (zi.le.0.00)  go  to  60 

c  if  (zl.le.0.00)  the  turning-point  has  been  reached. 

dtau(r)  =  2.0*sqrt(chl2)*om*x5*dh(r)/(xc*s(irt(zi)) 
if  (dtau(r)-l75.0)  25,26,24 


24 

ar(r)=l .00 

go  to  10 

25 

ar( r)=l .00- 

exp(-Jtau(r)  ) 

10 

con  t  1  nu  e 

go  to  61 

6  0 

p=r 

c 

p  =  pointer 

to  tlie  level  which  reflects  the  radio  wnvc  . 

c 

1’  C  r 

c  O'  i ,)  II  1  r*  ') ;)  t 

leal  ioptlis  of  direct  and  reflected  rays. 

6  1 

p  p  1  =  p  f  1 

Jr(rmpl )=0.00 
dtau(rmpl)=0.00 
j  t  »  rmp 1 - p 
do  30  j“l  ,  jt 
r  =  r  in  J)  1  -  j 
r  p  1 «  r+ 1 

30  dr(r)"‘dr(rpl)+dtau(rpl) 
c  dr(r)  •  optical  depth  of  level  r  (direct  ray). 

62  pp2“p+2 

do  27  r>«pp2,rm 

er(r)-dr(p) 

rml-r-l 

do  27  q=ppl,rml 
27  er(r)«er(r)+dtau(q) 

er(r)  -  optical  depth  of  level  r  (reflected  ray). 


c 


-  13  3  - 

er(ppl )=dr(p) 
tpuref(l)  =  ilr(p) 

c  tauref(i)  =  optical  depth  of  turnln'^-poiiit. 

c 

cccc  compute  contribution  from  each  level, 
do  32  r=ppl,rm 
b  r  (  r  )  = !)  .  0 

if  (d  r(  r  )-l  75 .0)  31,  2d,  2.3 

31  br(r)=br(r)Texp(-dr(r)) 

23  if  (cr(r)-175.())  36,32,32 

3  6  br(r)  =  br(r)3-exp(-er(r)) 

32  xr(r)  =  t(r)*.ir(r)*br(r) 

c  xr(r)  =  contribution  from  level  r  (k). 

c 

cccc  find  the  b  r  1  v,  t  .  temp.,  the  peak  of  t  lie  cont.  fn.  and  its  loc.it  ion, 
do  40  r=npl,rm 
cr(r)=xr(r)/dh(r) 

c  cr(r)  =  contribution  function  (k/km  at  level  r). 

if(crin(i)-cr(r))  33,40,4  0 
33crm(i)=cr(r) 

c  crm(i)  =  maximum  value  of  cr(r). 

r  c  rm  (  i  )  =  r 

c  rcrn(I)  =  location  of  the  peak  value  of  cr(r). 

4  0  bt(l)  =  ht(i)'t-xr(r) 

c  bt(l)  =  brinhtness  temperature  at  wavelenght  l(i)  (k). 

c 

cccc  normalize  the  contribution  and  tlie  contribution  function, 
c 

if  (  c  rm  (  1  )  .  oq  .  0 . 00  )  ,>o  to  72 

X  a  =  I  ,  0  /  c  r  m  (  I ) 
do  4  1  r  =  p p 1  , rm 
41  cr(r)=cr(r)*xa 

c  cr(r)  is  now  tiie  normalized  contribution  function. 

c 

72  xa-0.0 

do  69  r”ppl,rm 
if  (xa-xr(r))  63,69,69 
63  x.a“xr(r) 

6  9  c  o  n  t i nue 

c 

if  (xa.eq.O.OO)  go  to  73 
xa • 1 . 0 / xa 
do  6  7  r»pi)l,rm 
67  xr(r)»=xr(r)*xa 

c 

cccc  define  variables  for  output. 

7  3  a  r ( rnp 1 ) «0 . 00 

b  r ( rm  p 1 ) “ 1 . 00 
cr(rmpl  )*0.()0 
er(rmpl  )“0.0f) 
xr(rmpl  )»0.0() 
d  t  au  (  (>  )  ”0 . 0 
ar(p)-0.()0 
br ( p  )  *0 . 00 

c  r  (  p  )  “  0 . 0  0 


c 

c 

c  c  c  c 


/ 1) 
c 

c 

c  c  c  c 
4‘) 


-*  h 

c 

74 

c 

4  7 

4;i 

c 

c 

)') 

5  1 
c 

52 

5  1 

5  4 

5  5 

c 

c 

c  c  cc 


-  15  4  - 


1-- 1  (.  p  )  =  0 . 0  0 

x  r  (  p  )  =1)  .00 
L  lie  (  |)  )  =  U  .  0 

U  r  (  i  )  =  li  (  P  ) 

Ur(i)  =  help,  lit  of  reflection  level. 

com|)iite  theta  at  each  layer. 

y=»almpar*crd/r8un 

r  t  =  rra+-p  p  1 

do  7  0  j  =  ppl,rtn 

r=  r t - j 

x=y*dh(r)/sqrt(ro2(r)*(mu2(r)*ro2(r)-d2)) 
tho(r)  =  t!ic(rH  )+x 
r  o  n  t  i  n  11  e 

tlio(rjTpl)  was  def  tiled  above. 

L  h  e  i  11  f  (  i  )  =  L  he  (  p  p  1  ) 

deteriniiie  tlio  altitudes  of  the  half  power  points  of  cr(r). 
if  (crn(i).e<i.0.0  0)  qo  to  74 
j  =  rcrin(  i  ) 

i  ■'!  1  =  ,i  - 1 

do  4  4  r  =  n  ,  i "1 1 
s^  )  +  p-  1- r 

if  ( c  r ( s  )  -  0 . 5  )  4  5,47,46 

c  on  t  i  nil  e 

he  r2  (  I  )  =  (-  )<n  .0) 

he  r2  (  i  )  =-909 .0  means  tliat  the  model  does  not  go  deep  enoug.h. 

■’ll  to  50 

:icr2(  i  )  =  (-509.0) 
hcrl  (  l)  =  (-:UM.0) 

hcr-''(i)  =  -d.'j.'j.i)  means  there  is  no  contribution  at  this  wavelen  t  n . 

o  to  5  6 
!ic  r2  (  i  )  =  li(  s  ) 

■  o  to  5  0 
r  =  s  1 1 

hcr2(  i  )  =  ((0.5-cr(r))*(h(s)-h(r))/(cr(s)-cr(r)))  +  li(r) 
hcr2(i)  =  lower  half-power  point  of  the  cont.  fn. 

if  (j-rminl)  52,51,51 
(ic  r  1  (  i  )  =  (  -  9  '  9 . 0  ) 

lie  r  1  (  1  )  = -9  9  9  .  f)  means  that  the  model  does  not  go  high  cTiough. 

;;  o  to  5  6 
Jo  55  r=j,rmpl 
s  r 

if  ( c  r ( s  )  -  0 . 5  )  5  5,54  ,  53 

cont  ill u e 

hcrl ( i )=(-999.0) 
go  to  56 
licrl(l)=li(s) 
go  to  56 
r  =  s-  1 

hcrl  (  i  )  =  ((0.5-cr(r))*(h(s)-h(r))/{cr(s)-cr(r)))  +  li(r) 
hcrl(i)  =  upper  half-power  point  of  the  cont.  fn. 

output  contribution  function  for  wavelenp, th  l(i)  ccccccccccccccccc 


56  if  (In.nc.l)  Ko  to  64 

57  write(ll,616)(coin(J),j*=l,4()),rn,olf.i(l,b(>ta(),zota(),rl,alf.il, 

1  beta]  ,  z  e  t  a  1  ,  r  2  ,  a  1  f  a  2  ,  b  e  t  a  2  ,  z  e  t  a  2  ,  r  3  ,  .i  1  f  a  3  ,  !>  e  t  .t  3  ,  e  t  a  5  , 

2  "an  Tia  ,  a  l:npa  r 

uri te  (  1  1  ,  607  )  1  (  i  ) 

60/  format  (5x/16h  wavelcni;tti  =,10. 3, 56  (  c  m)  ,  64x  ,  2  5!iront  r  i  bu  t  i  on  f 
Iroin  level  r// 

1  1  3  )  li  r  1 1  cm  mo  *  *  2  d  t  a  o  a  r  dr  c  r 

2  b  r  t  h  e  t  a  0 . 0  0.2  0,4  0.6  l  .  K  1  . 

30  / 

4  1301,  (-)  (km)  (k/km)  (---)  (--)  (--)  (--)  (-- 

5)  (--)  (jee.) 

6  /) 
c 

do  34  r=p,rmpl 
kmax=xr(r)*40.0+l .001 
if  (  k  m  a  X  .  t  .  4  1  )  k  ma  x  =  4  1 

write  (  1  I  ,603  )  r,li(r),cr(r),i.iu2(r),dtau(r),ar(r),dr(r)  ,er(r), 

1  br(r),the(r),(ocr(k),k=l,kmax) 

60J  format  ( i 5 , f 3 . 0 , f 7 . 4 , f 7 . 4 , 5e 1 0 . 3  ,  f 7  .  2  ,  2 x  ,  4  1  i  1  ) 

if  (r.eT.rl.or.r.eq.r2.<>r.r.eq.r3.or.r.eq.r0)  write  (11,6  15) 

34  continue 
c 

cccc  find  the  heisht  at  which  tau  equals  unity. 

4  4  r  t • rm- p 

do  39  r»l,rt 
s“rm-r 

if  (dr(s)-l.O)  39,42,43 

39  continue 

htel(i)=(-999.0) 

c  htel(i)  =  -999.0  means  tau=l  not  reached  above  turninr' -point. 

1’, o  to  9  9  9 

42  htel(i)=h(s) 
jj  o  to  9  9  9 

43  r»s+l 

htel(i)=((1.0-dr(r))*(h(s)-h(r))/(dr(s)-dr(r)))+h(r) 
c  htel(i)  =  hei);ht  at  tau  =  l  at  wavelen.’tli  1(1)  (km). 

999  continue 
c 

cccc  wavelength  dependent  parameters  ccccccccccccccccccccccccccccccccc 
c 

lf(lp.eq.-l),’'o  to  312 

wrlte(ll  ,616)(com(j),j»l  ,40),r0,alfa0,beta0,zeta0,rl  ,alfal  , 

1  betal  ,zetal  ,r2,alfa2,bet.i2,zeta2,r3,alfa3,beta3,zeta3, 

2  gamma , a Impa r 
312  do  997  l-l,lr.i 

nu(i)»nu(i)*x9 
s»rcrm(  i) 
lnu-alop,l0(l(i)) 

1 f ( 1 p . e q . -  1 ) go  to  310 

997  write  (11,610)  l(i),nu(l),bt(l),lnu,h(s),htel(l),crm(l),hcr2(i), 

I  hcrl(l),hr(l),tauref(i),thelnf(l) 

format  (  5x  ,  I  4hwavc  lenp.t  h  =  ,  0  p  f  I  0 . 3  ,  2  6  h  cm  *  frequency 

l»,  f9.3,25h  ghz  *  brlgt.  temp.  =*,  fl0.(),7h  k  */ 


6  10 


-  136  - 


c 

iUO 

c 


(,  1  6 


c 


c 

20 


3  0  3 


c 

c 

c 

c 

c 


3 

4 


c 


3x,14l\lfi>>  unvlngt,  =,  fl0.3,2  6li  —  *  iic.  cfn.  pk. 


3-,  fO.l 

,2  3h 

km  *  h  t • 

t  a  u  “  1.0  = 

,  f  1  0  .  1  , 

7h 

k.n 

*/ 

4 

3x  , 

1 4  h  cfn. 

pk 

.  =  ,  1  p  e  1  0 . 3 

,  2  6  h 

k  /  km 

* 

lit  s 

.  cfn.  =  0.3 

3  =  ,')pf  9  .  1 

,4h 

=  =  >  ,  f  9 . 1  ,  3  h 

km 

,23x, Ih*/ 

b 

3x  , 

1  4  h  re  f  1  e  r  t  . 

ht 

.  =  ,  f  1  0 . 1 

,  26h 

km 

* 

t  :\  IJ 

at  r  e  f  1  e  i:  L  . 

7= , 1 pc9 . 2 

,2  3h 

_  * 

t  h . 

reflect.  * 

,  Op f 10.3, 

7h 

(1  G  . 

*/) 

if  (  I  c  oil  n  t  .  e  q  .  0  )  i’,  o  to  3 

w  r  i  t  o  (  6  ,  0  1  6  )  (  c  oin  (  i  )  ,  j  =  1  ,  A  0  )  ,  r  0  ,  a  1  f  a 0  ,  be  t  a  0  ,  z e  t  a  0  ,  r  1  ,  a  1 1  a  1  , 

1  be  t  a  I  ,  ze  t  a  1  ,  r  2  ,  a  i  f  a 2  ,  be  t  a 2  ,  ze  t  a 2  ,  r  3  ,  a  1  f  a  3  ,  be  t  a 3  ,  z e  t  a  3  , 

2  a  ;n  n  a  ,  a  i  ;n  p  a  r 

format  (  1  li  I  /  4  x  ,  4  Oa  2  /  /  3x  ,  6 h  r  0  =  ,  i  3  ,  1  1  h  ,  a  1  f  a  0  =  ,  f  b  .  2  ,  1  1  li  ,  lif- 
ItaO  =,f5.2,llli  ,  zetaO  ”,f5.2/  3x,6h  rl  =,i3, 

1  nil,  a  1  f  a  l  =  , 

'•  t  3. 2,  I  III  ,  lie  t  at  =,f3.2,ini  ,  zeCal  =  ,f3.2 

1/  3x,()li  r2  =,i3,llh  ,  alfa2  =,f5.2,llli  ,  beta2  =,fj.2, 

•  lib  ,  zfta2  =,f3.2/  3x,6h  r3  »  ,  i  3  ,  1  1  li  ,  alfa3  =,f3.2, 

31  111  ,  beta!  =■  ,  f  5  . 2  ,  1  1  li  ,  zeta3  =,f5.2/  3x,9li  ’.am  :ia 

3  fO.l,  13  li  ,  ai.a  par  ’',f6.3/) 

call  spoctr  (bt,bt<i,spec,l,afl«>,ir!i,}'tltl) 
a  f  1 1;=  .  true. 

wr i t  e ( 6 , o  0  3 ) r 0 , r 1  , r2 , r 3 , 1 p , a  1 f  aO  ,  a  1  f  a  1  ,  a  1 f  a2 , a  1 f  33  ,  be  t  aO  , 

Ibotal  ,bcta2,bcta3,zetaO,zetal  ,zeta2,zeta3,>’af'na 
f  o  riiui  t  (  '  1  i  npu  t  new  parameters.  to  enter  a  new  spectrum'/ 

I'  or  aiming;  parameter  enter  -1  for  rO'/ 

2  '  L  li  <■  o  1  il  p  a  r  a  rue  t  e  r  s  are:'  ,  3  I  3 , 4  f  6 . 2  /  4  9  x  ,  4  f  6  .  2  /  4  9  x  ,  4  f  6  .  2  ,  f  3 . 2  ) 
read  (9,*)  rO,rl,r2,r3,lp,alfaO,alfal  ,alfa2,alfa3,bocaO, 

1  be  t  a  1  ,  be  t  a2  ,  be  t  a3  ,  ze  t  aO  ,  ze  t  a  1  ,  ze  t  a 2  ,  ze  t  a 3  ,  ;;a mma 

tl  (l;).lt.-l)afli*  =  .false. 
if  (lo.lt.-l)lp=-l 
if  (  I  p  .  K  t  •  1  )  a  f  I  ;■  =  .  f  a  1  s  e  . 
i  I  (  1  P  .  a  t  .  1  )  1  p  =  1 
i f ( ze I aO . 1 e . 0 . 0 ) ze  taO* 1.0 
if(zctal  .le.0.'))zetal»1.0 
i  f  (zeta2  .  Ie.().0)zeta2=l  .0 
i  f  (  z  e  t  a  3  .  1  e  .  0 . 0  )  z  e  t  <a  3  =  1  .0 
i  f  (  r  0  .  e  (  ,  -  I  )  j;o  to  3  0  1 

Ip  is  a  switcli  with  toe  following  nieaiiin,'s: 

lp*0  means  (i  r  1  n  t  only  tlie  summary  table  (brightness  temperature, 
lp«l  means  print  also  the  normalized  contribution  functions. 
lp»  -1  means  only  print  out  the  spectrum. 

1  c  ou  tl  t  =  1  cou  n  t  1- 1 
r,o  to  7 


initialize  tlio  quiet  sun  for  this  a  in  par. 

do  4  1  =  1  ,  i  ni 

btn(i)“bt(l) 

btq(i)  »  quiet  sun  brightness  temperature  (k). 
go  to  20 


c 


-  n7  - 


9  9  > 

stop 

e  nd 

/Inc 

s  pck 

/  d  a  t  i 

) 

Zinc 

s  c  h  s  pe  c 

subroutine  s  p  e  c  t  r  (  h  t  ,  b  t  q  ,  s  p  e  c  ,  x  I  ,  a  f  1 1; ,  i  in  ,  t  i  t  I  e  ) 
i  II  t  e  .;e  r  t  i  1 1  e  (  1  5  ) 

real  *4  b t  (  2 1 ) , b t q ( 2 1 ) , s pe c ( 2 1 ) , x !( 2 1) 
real*4  qbt(2  )),  llm'la(2')),pspec(2  3) 
integer  cliar/O/ 
logical  a  f  I  g 
w  r  1 1  e  (  <1  ,  9  9  ) 

99  format('0  I  I md a ' , 5 x , ' r a t i o '  ) 
d  o  1  0  i  =  1  ,  I  m 

if(btq(i).eq.9.())btq(l)  =  ht(i) 
qbt(i)  =  bt(i)/htq(i)* lOO  ,0 
1  I  n  d  a  (  i  )  =  ,i  1  o  g  1 ')  (  x  1  (  i  )  ) 
uriLo(6,lii'))llnida(i),ql)t(i) 
l')f)  formatC'  ',2fl').4) 

1  ! )  continue 

char=cbar+l 

if  (  c  b  a  r  .  L  e  .  1  3  )  o  to  19 
a  f  1  g= . f  a  1 s  c . 
w  r  i  t  e  (  b  ,  1  0  1  ) 

101  for:nat('  cbango  plotter  paper') 

15  call  p  1  o  t  s ( 0 , 0 , 1 3 ) 
call  plot(l.i),1.0,-3) 
i  f  (  a  f  1  <’ )  >  o  to  50 
cha  r= I 

call  symhoKl  .5, 5. 0,0. 14, title, 0.0, GO) 
q  b t  (  i in+l  )  =  2 0 .0 
qbt  (  1  .^+2  )  =  20 .0 
Jo  11  i  =  1  ,  i in 
pspec(l)=spec( 1)*100.0 
1  1  lf(pspec(i).eq.0.0)pspec(l)*100. 
p  s  pe  c (  i  ra+  1  )  =  2  0 . 0 
p  s  pe  c  (  i  :t\+2  )  =  2  0 . 0 
call  s  c  a  1  e  (  I  I  111  d  a  ,  9  . 0  ,  i  m  ,  1  ) 

call  a  X  I  s  (  0 . 0 , 0  . !)  ,  1  0  h  I  og (  1  a.nda  )  ,  -  I  0 , 9 , 0 , 0 . 0 , 1 1  .nd  a  C  I  n  +  1  )  , 

11 ImdaC im+2) ) 

call  a xis(0.0,0.0,28h brightness  temi)erature  ratio, 2 3,5.0,9  0.0, 
1  pspec(  im+l  )  ,  pspec(  tr.i+2)  ) 
cal  I  line (I lmda,pspec, Im,  I  ,-l  ,0) 

50  call  llne(ll!iida,qbt,in),l,-l,char) 
call  plotf(().0,0.0) 
return 
e  nd 


138 


Sample  Input  and  Output  at  Computer  Terminal 


o 

c 


tn 

a 

IP 

X 

w 

X 


o 

CO 


o  o 


1 

CO 


ir. 

ts 

tu 


o  o  o 

00  GO  O 


c  o 


Ij' 


o  c  o: 
o  o  ♦ 


S3 

C-i 

uo 


X 

h-  c 

5J  so 


CO 

3 

u 

3 


2 

z> 


oo 


o 

-o 

■V 

r 

U ' 
IV 
I 

<2 


oo 


o  o  o  o 
iiO  o  c  o 

•  ♦  *  ♦ 

c-;  c-J  so  -ri 

M  II  i!  II 

o  ^  ri  ro 
<z  z.  z  < 

Lw  Li  iii  Ui 

^i  Ni  MM 


O  D  UT  O 

O'  oo  <r  03 

•  •  ♦  * 

o  o  o  o 

li  II  I)  II 

c  c  ;  r-o  o 
■zz  <:  -c  o 
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